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ABSTRACT 

The  problem  analysed  is  that  of  two-dimensional  wave  motion  in  a 
heterogeneous/  inviscid  fluid  confined  between  two  rigid  horizontal  planes  and 
subject  to  gravity  g.  It  is  assumed  that  a  fluid  of  constant  density  p+ 
lies  above  a  fluid  of  constant  density  p_  >  p+  >  0  and  that  the  system  is 
nondif fusive.  Progressing  solitary  waves/  viewed  in  a  moving  coordinate 
system  can  be  described  by  a  pair  (X/W)  where  the  constant  X  “  g/c2f  c 
being  the  wave  speed/  and  where  w(x/n)  +  n  is  the  height  at  a  horizontal 
position  x  of  the  streamline  which  has  height  n  at  x  =  ±».  It  is  shown 
that  among  the  nontrivial  solutions  of  a  quasilinear  elliptic  eigenvalue 
problem  for  (X,w)  is  an  unbounded  connected  set  in  R  x  (Hq  ft  C0,  1 ) . 

Various  properties  of  the  solution  are  shown  and  the  behavior  of  large 
amplitude  solutions  is  analysed  leading  to  the  alternative  that  internal  bores 
must  occur  or  streamlines  with  vertical  tangents  must  occur. 


AMS  (MOS)  Subject  Classifications:  35J60,  76B25,  76C10 

Key  Words:  Internal  wave/  Solitary  wave,  Internal  bore/  Bifurcation 
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SIGNIFICANCE  AND  EXPLANATION 


The  notion  of  a  solitary  wave  arose  in  the  last  century  and  provoked  some 
controversy  as  to  the  existence  of  such  a  phenomenon .  The  last  few  decades 
have  witnessed  an  upsurgence  of  scientific  work  on  these  waves  owing  to  the 
special  role  they  play  in  the  evolution  of  general  disturbances.  Solitary 
internal  waves  can  be  observed  in  the  laboratory  as  well  as  in  the  earth's 
oceans  and  estuaries  where  density  stratification  occurs,  due  to  changes  in 
temperature  or  salinity. 

The  problem  analysed  in  this  report  is  that  of  two-dimensional  wave 
motion  in  a  heterogeneous,  inviscid  fluid  confined  between  two  rigid 
horizontal  planes  and  subject  to  gravity  g.  It  is  assumed  that  a  fluid  of 
constant  density  p+  lies  above  a  fluid  of  constant  density  p_  >  p+  >  0  and 
that  the  system  is  nondiffusive.  Progressing  solitary  waves,  viewed  in  a 
moving  coordinate  system  can  be  described  by  a  pair  (X,w),  where  the 
constant  X  *  g/c2,  c  being  the  wave  speed,  and  where  w(x,n)  +  n  is  the 
height  at  a  horizontal  position  x  cf  the  streamline  which  has  height  n 
at  x  =  ±».  It  is  shown  that  among  the  nontrivial  solutions  of  a  quasil inear 
elliptic  eigenvalue  problem  for  (X,w)  is  an  unbounded  connected  set  in 
■  x  Hq  n  C0'1).  Various  properties  of  the  solution  are  shown  and  the  behavior 
of  large  amplitude  solutions  is  analysed  leading  to  the  alternative  that 
internal  bores  must  occur  or  streamlines  with  vertical  tangents  must  occur. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MFC,  and  not  with  the  authors  of  this  report. 


A  GLOBAL  THEORY  OF  INTERNAL  SOLITARY  HAVES  IN  TWO-FLUID  SYSTEMS 


C.  J.  Amick  and  R.  E.  L.  Turner 


1 .  Introduction 

The  study  of  single-crested  progressing  gravity  waves  was  initiated  over  a  century 
ago.  It  began  with  the  observations  by  Russell  [24)  of  what  he  termed  solitary  waves, 
which  progressed  without  change  of  form  over  a  considerable  distance  on  the  Glasgow- 
Edinburgh  Canal.  The  mathematical  analysis  of  this  wave  motion  on  the  surface  of  water, 
begun  in  the  nineteenth  century,  has  undergone  a  rapid  development  in  the  last  three 
decades,  due  to  the  scattering  theory  for  the  Xorteweg-de  Vries  equation,  which  models  the 
motion  of  long  waves  and  due  to  the  development  of  techniques  in  nonlinear  analysis 
allowing  for  the  analysis  of  finite  amplitude  motions. (cf.  [5),  [19]  and  references  given 
there ).^;  The  work  on  surface  waves  has  many  parallels  in  the  study  of  waves  in  fluids  with 
variable  density.  In  the  case  of  a  heterogeneous  fluid  with  a  free  upper  surface,  gravity 
waves  still  occur.  In  analogy  with  surface' waves  in  a  fluid  of  constant  density. (cf.  [23], 
[25],  [27] ) .  'What  is  distinctive  about  a  fluid  with  density  stratification,  however,  is 
the  presence  of  waves  which  are  predominantly  due  to  the  stratification  and  not  to  the 
free  surface.  These  waves,  called  internal  waves,  exist  in  a  heterogeneous  fluid  even 
when  it  is  confined  between  horizontal  boundaries,  a  configuration  which  precludes  gravity 
waves  in  a  fluid  of  constant  density. 

For  surveys  of  earlier  work  on  permanent  waves  in  stratified  fluids  and  for  more 
complete  references  than  given  here,  we  refer  the  reader  to  the  articles  by  Benjamin  [8], 

Bona,  Bose,  and  Turner  [10],  and  Turner  [26]. 

.  >■ ■  _ 

Our  concern  in  this  paper  is ,with  progressing  solitary  gravity  waves  in  a  system 
consisting  of  two  fluids  of  differing  densities  confined  in  a  channel  of  unit  depth  and 
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infinite  horizontal  extent.  The  undisturbed  state  consists  of  a  fluid  of  constant  density 
p_  having  depth  h  in  contact  with  the  lower,  planar  boundary  and  a  second  fluid  of 
constant  positive  density  p+  <  p_,  lying  above  the  first  and  having  depth  1  -  h,  so  as 
to  fill  the  channel.  The  motion  studied  is  assumed  to  be  two-dimensional,  inviscid, 
incompressible,  and  nondiffusive.  We  shall  always  remove  the  time  variable  from 
consideration  by  taking  a  moving  coordinate  system  in  which  the  flow  is  steady  and 
symmetric  about  the  vertical  axis. 

The  problem  of  finding  nontrivial  wave  motion  can  be  posed  in  terms  of  a  pair  of 
harmonic  stream  functions,  1>+,  for  the  upper  fluid  and  t~,  for  the  lower  fluid.  At 
the  interface  between  the  two  fluids,  a  free  boundary  in  this  setting,  one  must  impose  the 
continuity  of  pressure  and  of  normal  velocity  components.  A  direct  analysis  of  the 
problem  in  this  setting  might  be  attempted  along  the  lines  of  the  work  on  two  fluid  jets 
by  Alt,  Caffarelli,  and  Friedman  ((1],  [2]).  However,  they  make  essential  use  of  having  a 
minimizer  for  an  appropriate  functional  and  must  do  a  delicate  analysis  to  show  that  the 
interface  separating  the  two  fluids  is  of  class  C1.  In  the  present  context  the  analogous 
minima  would  correspond  to  trivial  parallel  flows  <cf.  [10],  [26])  and  the  nontrivial 
flows,  to  what  are  undoubtedly  saddle  points.  By  taking  a  different  approach  we  obtain  a 
global  picture  of  the  solitary  wave  solutions  and  show  that  for  each  wave  the  interface 
between  the  two  fluids  is  analytic. 

The  approach  taken  in  this  paper  is  to  approximate  the  discontinuous  density  function 
by  a  sequence  of  smooth  density  profiles  and  to  obtain  sufficiently  good  estimates  in  the 
smooth  case  to  allow  appropriate  limits  to  be  taken.  In  this  way  we  also  show  that 
discontinuous  densities  and  smooth,  rapidly  varying  densities  give  rise  to  motions  which 
are  close  to  each  other.  Using  methods  of  global  bifurcation,  the  first  author  in  [4] 
gave  a  global  theory  for  solitary  waves  in  the  case  of  an  arbitrary  smooth  density.  The 
solutions  are  given  by  stream  functions  which  are  smooth  everywhere  in  the  channel.  This 
approach  is  not  limited  by  the  particular  shape  of  the  density  profile,  as  was  the 


variational  approach  uaad  in  [10].  However,  a  limiting  process  to  obtain  solutions  for 
the  density  having  a  jump,  from  p_  to  p+,  would  still  entail  a  confrontation  with  an 
unknown  free  boundary. 

An  alternative  is  to  use  semi-Lagrangian  coordinates.  To  describe  these  coordinates 
we  first  note  that  the  stream  function  referred  to  in  the  last  paragraph  differs  from  the 
usual  one  in  having  an  intervening  density  factor  (see  formula  2.14).  Using  accepted 
terminology  we  call  it  a  pseudo-stream-function  and  denote  it  by  ^  «  4>(x,y),  where  x 
is  a  horizontal  coordinate  and  y  is  a  vertical  coordinate  chosen  so  that  gravity  acts  in 
the  negative  y  direction.  In  the  semi-La grangian  formulation  one  uses  x  and  as 
independent  variables  and  y  as  a  dependent  variable.  In  this  formulation  the  equation 
to  be  considered  is  highly  singular.  However,  what  would  be  free  boundaries  in  the 
spatial  domain  correspond  to  coordinate  lines,  4>  *  constant,  in  the  new  variables.  This 
formulation  was  used  by  Ter-krikorov  [25]  in  combination  with  fixed  point  methods  to  show 
the  existence  of  small  amplitude  periodic  and  solitary  waves  in  the  case  of  a  fixed  or 
free  boundary  at  the  upper  surface  of  a  fluid  with  smoothly  varying  density.  Using  this 
formulation,  but  applying  a  variational  method,  the  second  author  in  [26]  and  [27]  studied 
the  two  fluid  system  under  consideration  here  as  well  as  a  multiple  fluid  system  with  a 
free  upper  surface.  The  existence  of  periodic  and  solitary  waves  was  shown,  the  latter 
waves  being  obtained  as  the  limit  of  periodic  ones  with  ever  increasing  periods.  In  the 
present  paper,  to  obtain  large  amplitude  waves  in  the  case  of  a  discontinuous  density  we 
work  principally  with  the  semi-Lagrangian  formulation,  but  as  a  starting  point  use  the 
solutions  provided  by  [4].  It  will  be  shown  that  there  is  a  global  theory  for  a 
discontinuous  density  analogous  to  that  derived  in  [4]  for  a  smooth  density. 

In  Section  2.  1  the  passage  from  the  postulated  physical  model  to  the  relevant 
equations  is  made.  The  equations  are  of  two  types.  The  first  is  a  semi-linear  elliptic 
equation  for  a  pseudo-stream-function  ^(x,y).  The  second  is  a  singular,  quasi-linear 
elliptic  equation  which  governs  the  deviation  of  streamlines  from  parallel  flow  when  semi- 
Lagrangian  coordinates  are  used.  In  both  equations  there  enters  an  unknown  parameter 


t 


which  describes  the  speed  c  of  a  progressive  wave  as  seen  in  "laboratory"  coordinates. 

2 

In  the  second  formulation  it  appears  explicitly  as  an  eigenvalue  parameter  X  =■  g/c 
where  g  is  the  gravity  constant. 

In  Section  2.2  the  main  results  of  the  paper  are  stated  in  Theorem  2.1.  It  describes 

a  branch  S  of  solution  pairs  (X,w)  where  w  -  w(x^,x2)  is  the  deviation  at  the 

horizontal  position  Xj  of  the  streamline  which  has  height  x2  at  x^  “  ±<*.  If  T 

denotes  the  doauin  of  the  pairs  (x1.x2) ,  the  branch  is  unbounded  and  connected  in 
i  o.l  - 

ft  x  (Hq(T)  n  C  (T) ) .  It  emanates  from  a  pair  (1^,0)  where  X^  is  a  particular 

1/2 

function  of  p+,  p_,  and  h.  The  speed  c^  -  ( g/X^)  ie  a  critical  speed  for  the 
stratified  configuration,  the  analogue  of  the  speed  /gh  for  long  waves  on  the  surface  of 
water  of  depth  h.  For  all  nontrivial  solutions  in  S,  0  <  X  <  X^  and  thus  the 
associated  speeds  are  supercritical,  a  result  consistent  with  all  analyses  of  solitary 
wave  phenomena. 

The  sign  of  another  parameter  e,  a  function  of  p+,  p_,  and  h,  predicts  the  sign 
of  streamline  displacements.  If  e  >  0,  the  undisturbed  state  will  allow  waves  of 
elevation,  while  for  e  <  0  it  allows  waves  of  depression.  The  phenomenon  described 
here,  as  distinct  from  that  associated  with  simple  eigenvalue  bifurcation,  is  one-sided. 

If  e  >  0  there  are  elevation  waves  but  no  depression  waves  near  the  bifurcation  point 
(X^.O).  Finally,  all  waves  are  symmetric  about  their  crests  and  decay  monotonically  and 
exponentially  to  zero  as  x  ■*  ±«*. 

The  core  of  the  proof  of  Theorem  2.1  is  carried  out  in  the  course  of  Sections  3 
through  6.  We  begin  with  a  sequence  pr>  n  -  1,2,3,...,  of  smooth  density  profiles  in 
the  undisturbed  channel  configuration  which,  as  n  ♦  •,  converge  to  the  given 
discontinuous  density  function.  It  is  shown  in  Section  3  that  for  each  n  a  subset  of 
the  peeudo-stream-functions  from  (4]  provide  an  unbounded,  closed  connected  set  of 
solutions  (X,w)  in  R  x  Hq(T)  with  which  to  work  in  the  new  coordinates.  In  Section  4 
a  general  scheme  for  the  limiting  process  is  given.  Then,  after  imposing  restrictions  on 
the  gradient  of  w,  regularity  estimates  for  solutions  are  derived,  the  estimates  being 
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independent  of  n.  These  estimates  are  used  in  Section  5  to  obtain  connected  sets  of 
solutions  in  R  x  Hq(T)  for  the  case  of  a  discontinuous  density.  In  Section  6  the 
solution  sets  are  examined  in  the  stronger  topology  of  R  x  (Hq(T)  nc0,1(T)).  Here  it  is 
shown  that  connectedness  is  preserved.  The  next  step  concerns  the  restrictions  imposed  on 
the  gradient  of  w  in  Section  4.  It  must  be  shown  that  for  any  solution  (A,w)  under 
consideration,  the  gradient  of  w  does  not  take  values  in  a  range  that  would  make  the 
underlying  elliptic  equation  singular.  This  is  accomplished  by  returning  to  tKs  original 
formulation  involving  the  pseudo-stream-functions  and  obtaining  bounds  on  their 

gradients;  that  is,  bounds  for  the  velocities  in  the  flow  domain.  These  bounds  translate 
into  the  desired  estimates  for  the  gradient  of  w. 

In  Section  7  results  constituting  the  proof  of  Theorem  2.1  are  collected.  Many 
properties  of  solutions  are  obtained  in  the  course  of  proving  the  existence  of  the  set  5 . 
To  these  are  added  the  proofs  that  the  interface  is  analytic  and  that  depression 
(elevation)  waves  are  absent  near  the  bifurcation  point  (Ad,0)  if  e  >  0  (e  <  0). 

Numerical  studies  have  been  made  of  periodic,  interfacial,  gravity  waves  for  two 
unbounded  fluids  of  differing  densities  (cf.  [17),  [28]).  It  is  found  that  along  a  branch 
of  waves  of  fixed  wavelength,  the  separating  streamline  steepens,  eventually  manifesting  a 
vertical  tangent  and,  past  that,  an  overhanging  region  in  which  heavier  fluid  lies  above 
the  lighter  one.  From  these  studies  one  might  conjecture  that,  along  the  branch  of 
solutions  S  found  in  this  paper,  w  becomes  large  in  C0,1.  In  fact,  the  introduction 
of  the  C0'1  topology  apart  from  having  a  crucial  role  in  the  estimates  of  Section  4,  was 
suggested  by  the  numerical  work.  Section  7  concerns  the  behavior  of  large  amplitude 
solutions  on  the  branch  S.  The  fact  of  having  a  channel  of  finite  depth  and  solitary 
rather  than  periodic  waves  may  change  the  character  of  large  amplitude  waves.  He  pursue 
the  implications  of  assuming  that  wave  profiles  do  not  steepen  to  the  point  of  having 
vertical  profiles.  A  first  consequence  is  that  there  must  exist  a  solution  (A,w)  of  the 
flow  equations  for  which  w  is  not  in  L  (T).  Rather,  at  x  *  £■  it  asymptotically 
approaches  a  nontrivial  parallel  flow,  a  "conjugate”  flow  in  the  terminology  of  Benjamin 


[9] .  A  further  consequence  is  the  existence  of  an  "internal  bore"  of  predictable  size  am 


speed.  As  yet,  no  contradiction  has  arisen  from  this  train  of  arguments  and  so  the  large 
amplitude  behavior  remains  an  open  question.  In  a  project  in  progress  the  second  author 
and  J.-M.  Vanden-Broeck  are  carrying  out  a  numerical  study  of  the  solitary  wave  patterns 
shown  to  exist  here  and  hope  to  shed  light  on  their  behavior. 
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2.  Formulation  of  the  Flow  Problem  and  Results 
2.  1.  The  Flow  Equations 

Here  we  describe  the  physical  model  under  consideration  and  the  passage  to  a  boundary 
value  problem  for  a  partial  differential  equation.  The  boundary-value  problem,  in  turn, 
is  formulated  in  essentially  two  ways  and  the  interplay  between  these  is  used  in  our 
analysis.  Consider  an  incompressible,  inviscid  fluid  acted  on  by  gravity  and  restrict 
attention  to  a  two-dimensional  flow  confined  to  and  filling  a  region 

S  -  {(x,y)|-  <  x  <  ,  -h  <  y  <  1  -  h}.  The  acceleration  of  gravity  has  magnitude  g  in 
the  negative  y  direction.  Further  assumed  characteristics  of  the  fluid  which  make 
propagation  of  permanent  waves  possible  and  their  study  tractable  are  that  it  is 
heterogeneous  and  non-diffusive.  To  explain  the  last  term  it  is  worthwhile  t^>  first 
describe  a  diffusive  system.  Consider  a  mixture  of  two  ompcnents,  say  water  and  a 
dissolved  salt.  In  a  mixture  the  molecules  can  trade  places  without  having  any  net  mean- 
molar  velocity.  However,  if  the  molecules  have  different  masses,  there  can  be  a  net 
movement  of  mass.  Let  q  denote  the  mean-molar  velocity  and  Q,  the  mean-mass  velocity 
(cf.  [11],  Chapter  1).  Let  p  denote  the  density  of  the  fluid  and  suppose,  to  begin 
with,  that  p  is  a  smooth  function  of  space  and  time.  Subsequently  we  shall  formulate 
the  problem  for  a  discontinuous  density.  If  v  denotes  the  diffusion  coefficient 
(assumed  the  same  for  the  two  types  of  molecules)  then  the  Flck  diffusion  law  (cf.  [7]*> 
[11]  i  p.  23?  [14]  )  gives 


The  Euler  momentum  equation  is 

-a 

p(-|&  +  (Q  «  V )Q)  »  -Vp  -  pgk  (2.2) 

where  p  is  the  pressure  and  H?  Lb  a  unit  vector  in  the  y  direction.  Conservation  of 
mass  takes  the  form 

*We  wish  to  thank  the  authors  for  a  discussion  of  the  diffusive  model. 
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||  +  div(pQ)  =  0  (2.3) 

and  the  volume  incompressibility,  the  form 

div  'cf  ■  0  .  (2.4) 

A  set  of  equations  for  p  and  "q  is  obtained  by  using  (2.1)  in  (2.2)  and  (2.3).  We  will 
seek  a  wave  motion  which  has  a  permanent  form  when  viewed  in  a  coordinate  system  moving  in 
the  x  direction  with  constant  speed  and  in  such  a  coordinate  system  the  equations  will 
have  no  explicit  time  dependence.  To  avoid  introducing  new  notation  let  us  assume  for  the 
remainder  of  the  paper  that  all  quantities  are  referred  to  the  moving  system,  but  retain 
the  notation  already  introduced.  With  no  time  derivatives  present,  the  equations  for  p 
and  cf  reduce  to 

-A  -X  -A 

p(q  •  V)q  +  0(v)  =  -Vp  -  pgk  ,  (2.5) 

A 

q  •  Vp  -  vAp  «  0  ,  (2.6) 

div  f =  0  (2.7) 

where  O(v)  refers  to  terms  of  order  v  for  v  small.  Time  scales  for  diffusion  (e.g. 
in  stratified  salt  solutions)  are  on  the  order  of  hours  while  a  wave  of  the  type  sought 
here  will  pass  in  a  matter  of  seconds.  Thus,  diffusion  plays  an  insignificant  role  in  the 
wave  propagation  phenomenon.  By  a  nondiffusive  fluid  we  merely  mean  a  fluid  governed  by 
(2.5)-(2.7)  with  v  «  0.  The  resulting  equations  are  the  standard  ones  used  for 
nondiffusive  stratified  fluids,  but  the  condition  "q  •  Vp  «  0,  attributed  to 
nondiffusivity,  is  hereby  rationalized.  Otherwise,  it  would  seem  to  arise  from  using  a 
single  solenoidal  field  <f  together  with  mass  conservation  and  not  have  any  clear 
connection  with  diffusion  (or  lack  of  it). 

The  equations  which  result  from  setting  v  =  0  in  (2.5)— (2.7)  and  which  are  the 


basis  for  our  analysis  are : 


p(q  •  V)q  =  -Vp  -  pgk  , 
q'  •  Vp  =  0  , 


div  q  =  0 


(2. 10) 


For  use  in  the  sequel  we  let  q  *  (U,V)  where  U  and  V  denote  the  horizontal  and 


i] 

V 


\ 

N 


4 


vertical  components  of  velocity,  respectively.  It  will  be  assumed  that  in  the  original 


"laboratory"  coordinates  the  fluid  is  at  rest  at  x  =  and  there  has  a  positive  density 


P*(y)  which  is  nonincreasing  in  y  and  satisfies  p^t-h)  >  P„(1  -  h).  Thus  if  a  wave 


under  consideration  is  moving  from  right  to  left  in  the  "laboratory"  with  speed  c  >  0,  a 


corresponding  boundary-value  problem  in  the  new  coordinates  is  to  find  q(x,y). 


p(x,y),  p(x,y),  and  c  satisfying  (2.8)-(2.10)  and 


lim  p(x,y)  -  p  (y)  , 

l*l~ 


(2.11) 


lim  <U,V)  -  (c,0)  , 


(2.12) 


V  -  0  on  3S 


(2.13) 


where  3S  denotes  the  boundary  of  S.  Note  that  'q(x,y)  »  (c,0),  p(x,y)  =  -  /  pa(s) gde , 

0 


and  p(x,y)  »  PM(y)  is  a  solution  of  (2.8)-(2.13)  for  any  real  c  and  will  be  referred 


to  as  a  trivial  solution.  We  shall  find  nontrivial  solutions  of  the  boundary  value 


problem  ( 2.8)-( 2. 13 ) ,  properly  interpreted,  for  a  piecewise  constant  density  p^fy) , 


given  as  p  (y)  and  shall  do  this  by  way  of  smooth  approximating  densities.  In  the 


treatment  of  smooth  or  discontinuous  densities  it  will  always  be  assumed  that  each 


streamline  (i.e.,  integral  curve  for  q>  is  a  simple  Jordan  arc  connecting  x  «  -»  to 


x  «  +*■>.  In  particular,  there  are  to  be  no  "internal  eddies",  that  is,  no  relatively 


compact  regions  in  S  bounded  by  parts  of  streamlines. 


From  (2.9)  and  (2.10)  it  follows  that  q  admits  a  stream  function,  as  does  q 


multiplied  by  any  function  of  p,  since  p  is  locally  constant  on  streamlines  for  q. 


In  particular,  there  is  a  "pseudo-stream-function"  i|i(x,y)  for  which 

1/2  -»> 

a  ’ '  *  ~  —  /  ,1.  _.l.  \ 


q  -  ( ^  ' 


(2.  14) 


the  subscripts  denoting  partial  derivatives.  The  total  head  pressure  H  is  defined  by 


H(x,y)  -  p(x,y)  +  4-  p (x,y)  |q( x,y )  | 2  +p(x,y)gy 


(2. 15) 


and  from  (2.8)  it  follows  that 


U-  0V(VX  -  d  ,  ♦  PX(J  |q|2  -  gy> 


(2. 16) 


M  -  pU(Uy  -  Vj£)  +  p y(l  |q|2  +  gy)  .  (2.17) 

Since  ^  •  Vp  =  0  the  Bernoulli  condition 

q  •  VH  •>  0  (2. 18) 

follows.  Hius,  H  and  p  may  be  considered  as  functions  of  the  single  variable  41  and 
we  shall  write  H(x,y)  or  H(i|i)  and  p(x,y)  or  p ( 4*)  when  no  confusion  is  possible. 
Since  all  streamlines  go  to  infinity  where  p  and  q1  are  known  and  where  p  can  be 
assumed  hydrostatic: 


y 

p(y)  »  _  f  p^fsjgds  ,  (2.19) 

0 

the  functions  H(<|>)  and  p(i|i),  can  be  computed  explicitly  when  pm  and  c  are  given 
(cf.  (2.21),  (2.22)).  The  explicit  examples  in  [10],  section  6,  illustrate  the 
computation. 

Define 


V  _ 

T<y)  =  c  /  /p.(s)  ds  ,  (2.20) 

0 

the  pseudo-stream-function  "at  infinity".  Let  Y(t|i)  denote  the  function  inverse  to  T 
so  that  Y(¥(y))  =  y  for  -h  <  y  <  1  -  h.  Then  all  along  the  streamline  with  value  <i, 
the  density  has  the  value 

p(i|i)  -  P„<Y<*))  (2.21) 

and  for  a  flow  with  a  pseudo-stream-function  <i(x,y )  the  density  at  (x,y)  is 
P(x,y))  =  P„(Y(  <|i(x,y) ) ) .  That  is,  the  streamline  through  (x,y)  must  be  followed  to  «• 
to  ascertain  its  associated  density.  Similarly,  from  the  data  at  infinity  we  have 

2 

H(iji)  -  p(Y(«))  +  p(«){|-+  gY(4>)}  (2.22) 


and  each  of  (2.21)  and  (2.22)  is  defined  for  'f(-h)  <  *  <  Y(1  -  h) 


to  obtain  a  partial  differential  equation  for  suppose  U  f  0  (as  we  shall). 

Using  the  chain  rule  to  express  the  derivatives  of  H  and  p  in  (2. 17)  one  obtains 


S-  ^  ‘°v  -  V  <T  l£|2  ♦  w> 


d*  '2 


after  cancelling  a  factor  /p  U  «  If  U^  -  Vx  is  evaluated  using  (2.14),  the  last 


A* +  »  H  ■  If 


equation  becomes  the  semilinear  elliptic  equation 

(2.23) 

which  is  Yih's  version  of  Long's  equation  (cf.  [10]  for  references).  The  associated 
boundary  conditions  are 

4i(x,-h)  -  »<-h),  i(i(x,  1  -  h)  -  f(1  -  h >i  x  e  R  (2.24) 

and 

11m  *(x,y)  -  ?(y)  .  (2.25) 

|x|-*~ 


The  formulation  at  this  stage,  given  Pa(y),  is  to  find  c  and  a  function  f(x,y) 
satisfying  (2.23)-(2.25).  If  p^  is  smooth,  one  seeks  a  smooth  i|i.  If  pm  Is  the 
piecewise  constant  function 

p_  -h  <  y  <  0 


Pd<y>  -  i 


(2.26) 


p+  0  <  y  <  1  -  h 


then  the  meaning  of  (2.23)  must  be  clarified.  The  dividing  streamline  between  the 
portions  of  fluid  with  different  densities  is  that  on  which  ■  0.  Suppose  it  is  the 
graph  of  a  function  T  -  { (x,y  (x) )  |  -*•  <  x  <  •}.  In  each  of  the  regions 

S"  ■  {(x,y)|-h  <  y  <  y(x)}  and  S+  ■  {(x,y)|y(x)  <  y  <  1  -  h),  p  is  constant,  taking  the 
values  P  and  p+,  respectively.  From  (2.19)  and  (2.22)  H  is  constant  in  each  region 
and  by  calculating  its  values  at  x  »  *  and  y  »  0,  on  either  side  of  T,  one  sees  the 

values  to  be  H”  •  p_c2/2  and  H+  -  p+c2/2,  respectively.  Equation  (2.23)  requires  ♦ 

± 


to  be  harmonic  in  S  ,  as  expected,  ks 
continuous .  Thus 

♦  «  0  on  r 


is  a  measure  of  flux  we  require  it  to  be 


(2.27) 
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and  since  it  should  also  satisfy  (2.24) ,  by  the  maximum  principle  it  will  be  positive  in 
S+  and  negative  in  S”.  Let  denote  the  restrictions  to  S*,  respectively.  An 

additional  requirement  on  the  free  boundary  r  is  that  the  pressure 

P  *  H  -  j  |V#|2  -  pgy 

be  continuous.  Thus 

2 

\  {|?*+|2  -  I?*  |2}  -  <P  +  -  P_>{| - 9y>  on  r  (2.28) 

using  the  values  H*  from  above.  To  summarize,  the  boundary  value  problem  for  the  given 
discontinuous  density  (2.26)  is  to  find  T,  c,  and  f(x,y)  such  that  is  harmonic  in 
S*  and  satisfies  t2.24),  (2.25),  (2.27),  and  (2.28).  There  will  be  no  problem 
interpreting  the  conditions  as  it  will  be  seen  that  T  is  smooth  as  are  the  extensions  of 
to  the  closures  S*. 

Next  we  give  alternate  formulations  of  the  boundary-value  problems  corresponding  to 
the  smooth  and  discontinuous  densities,  respectively.  He  shall  work  only  with  flows  in 
which  no  reversal  occurs;  that  is,  for  which  U(x,y)  >0  in  S.  From  (2.14)  this 
corresponds  to  >  0  and  for  such  functions  one  can  solve  for  y  as  a  function  of  the 
spatial  coordinate  x  and  the  "material  coordinate”  <|».  The  utility  of  this  semi- 
Lagrangian  description  is  that  for  the  density  (2.26)  the  unknown  interface  function  y 
>is  merely  the  unknown  function  y(x,^)  evaluated  on  ^  =  0;  i.e. 

y(x)  -  y(x,0)  .  (2.29) 

A  disadvantage  of  using  y  as  a  dependent  variable  is  that  the  semi-linear  equation 
(2.23)  is  replaced  by  a  singular  quasi-linear  equation  for  y(x,<i).  Since 
y(x,f(x,y))  »  y  one  has  the  relations 

yx  +  Vx  *  ° 

V  (2.30) 

y .  •  *  -  1 

♦  Ty 

from  which  one  derives  the  equation 


1 1  s  «.  j"  vm  u  .■  sv  ww 


Bi 


n 


$ 


K 


in  It  «  (T(-h),T(1  -h) )  from  (2.23).  The  associated  auxiliary  conditions  ara 

y(x,T(-h))  -  -h,  y(x,T( 1  -  h) )  -  1  -  h,  x  e  R  (2.32) 

and 

lia  y(x,t)  -  Y«i)  .  (2.33) 

M** 

Ona  final  changa  of  variables  will  bring  the  aquation  to  tha  fora  used  by  the  aacond 
author  in  [26] .  Let  x1  ■  x,  Xj  »  Y($)  and 

wfx^Xj)  “  ytx^KXj))  -  x2  (2.34) 

ao  that  w  represents  tha  deviation  at  a  horizontal  position  x^  of  the  straaaline  which 
has  height  Xj  at  «.  it  should  be  emphasized  that  x2  is  a  rescaled  straaa 
coordinate.  In  describing  an  equation  for  w  we  set 


P,  *  P2 


(2.35) 


9f 


^Pi'P*5  *  2(1~  p2)  ' 

2 l  T  -  m  x  (-h, 1  -  h) >  and  use  the  summation  convention,  A 


f^  ”  i  i  -  1,2)  X  ■  g/e 
calculation  produces  an  eigenvalue  problem  for  the  pair  (X,w)> 


-  •j~-  Pmlx2)ti{V*)  -  -Xp^(x2)w  in  T  , 

(2.36) 

w(x,,-h)  ■  w(x,,1  -  h)  »  0,  x,  e  R  , 

(2.37) 

lim  w(x.,x2)  -  0  . 

(2.38) 

|x|*- 

The  case  for  the  discontinuous  density  (2.26)  requires  a  separate  formulation.  While  T 
and  8  coincide,  the  domain  T+  »  {(x^Xj)  e  t|x2  >  0}  is  what  corresponds  to  S+.  One 
defines  T“  similarly.  By  a  solution  of  (2.36)-( 2.38)  for  the  density  (2.26)  vr 
understand  a  function  w  e  C°' 1 (T)  n  C 1 ( T* )  satisfying  (2.37),  (2.38),  and  the  weak 
equation 


/  -  -X  /  p’»t  -  X(p_  -  p+>  /  w(x1#0)f (x1 ,0)dx1 

Til  -m 


(2.39) 
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for  all  f  e  Cg(T).  To  keep  both  (2.36)  and  (2.39)  in  aind  we  shall  aiaply  refer  to  a 
solution  of  (2. 36) -(2. 39) .  Throughout  the  paper  standard  notation  for  function  spaces 
will  be  used  (cf.  [13]).  He  shall  write  C  for  a  space  C°  of  continuous  functions. 


2.2.  Results 

Theorem  2. 1 .  Let  p^  be  the  discontinuous  density  given  in  (2.26),  let 


i  P_  P* 

-  _ . 1  .  , 

(p  -  p  )  '•h  1  -  hJ 


(2.40) 


and  let 


2  2  * 
hx  (1  -  h)* 


(2.41) 


Then  if  e  >  0, 

(a)  There  exists  an  unbounded,  closed,  connected  set  $C  R  x  (Hg(T)  n  C0,1(T))  of 
solutions  (X,w)  of  (2. 36) -(2. 39)  with  (1^,0)  t  S. 

(b)  There  is  a  positive  constant  X  -  X(p  ,p»  such  that  if  (X,w)  e  S\{ (X.,0)} , 

—  —  •  ▼  a 

9  v 

then  £  <  X  <  X^,  w  >  0  on  T,  wi-x^Xj)  -  w(x^,x2)  on  T,  and  <  0  on 

(0,«)  x  ( -h, 1  -  h). 

(c)  The  function  w  has  real  analytic  extensions  to  T*"  and  T~,  e  C°(T),  and 

8x1 


Pd<X2,{f2(7w)  "  Xw}  C  C  <T) 


(2.42) 


for  some  HSlder  exponent  a  >  aQ  >  0,  uniformly  for  w  in  bounded  sets  of 
C0,1(T). 

<d)  There  exists  a  constant  K  ••  K(p_,p+,h)  such  that  for  all  (X,w)  e  S, 

0  <  - - - -  <  K  in  T 

1  +  w 

X2 


<  K  in  T 


wm$ 


*  1 


(a)  For  Mch  (X,w)  e  S\{(Xd,0)} 


|w(x1fx2)|  +  |Vw{x1#x2>|  <  C,*  ^ 

where  C 1  and  C2  ara  positive  constants  depending  on  (X.w) . 

(f)  If  | X  -  X  |  +  |w|  .  +  | w (  .  is  sufficiently  small,  the  only  solution 

H0(*)  c  ' 

of  (2.36)-(2.39)  for  which  w  is  even  in  Xj  and  for  which  both  w  and 
-3w/3x  have  signs  opposite  to  that  of  e  for  >  0  is  the  trivial 
solution* 


If  e  <  0,  the  results  (a)-(f)  hold  with  the  change  w  <  0  in  T  and 
for  xj  >  0  in  part  b). 


We  shall  restrict  attention  to  proving  parts  (a)-(f)  for  e  >  0  since  the  case 
e  <  0  is  analogous* 

The  implications  of  the  assertions  in  the  theorem  for  the  fluid  flow  are  not  all 
transparent  and  soma  discussion  is  in  order.  The  value  X^  at  which  bifurcation  takes 
place  and  the  sign  of  e  as  an  indicator  of  the  direction  of  streamline  displacements 
both  arise  in  the  small  amplitude  calculations  of  bong  [16]  and  Benjamin  [B]  as  well  as  in 
the  variational  approach  of  the  second  author  in  [26] .  In  the  passage  from  the  physical 
flow  model  to  the  problem  (2.36)-(2.39)  we  did  not  address  the  question  of  which  solutions 
(X.w)  give  rise  to  physically  reasonable  flows.  In  the  approaches  using  the  equation 
(1.23)  for  f(x,y)  further  restrictions  Intervene  (cf.  [4],  [10]).  Here,  however,  each 
pair  (X.w)  €  S  corresponds  to  a  physical  flow.  The  main  point  is  to  verify  that  all 
streamlines  connect  x  «  -*•  to  x  "  +"  or,  equivalently,  that  -h  <  x2  +  w(x1,x2)  <  1  -  h 
for  -h  <  x2  <  1  -  h.  Since  w  •  0  for  x2  “  -h  and  x2  *  1  -  h  and  since  1  +  wx^  >  0 
by  part  d),  it  follows  that  x2  +  w  lies  in  the  desired  range.  The  bounds  in  part  d) 
are,  in  reality,  bounds  on  velocities,  for  the  use  of  the  relations  (2.30)  together  with 
(2.14)  and  (2.34)  yield 


-v»v *.*.v  v  v  ■  • . 


The  velocity  c  satisfies 


(g/X.)1/2  <  c  <  <g/X) 


1/2 


(2.44) 


using  part  b)  and  thus  U  and  V  have  bounds  depending  only  on  p_,  p  +,  and  h.  The 

velocity  (g/X.)1^2  is  the  so-called  "critical  velocity"  associated  with  the  fluid  system 
a 

and  (2.44)  shows  that  all  solitary  waves  found  here  have  supercritical  velocities.  The 

remaining  parts  of  b)  describing  the  shape  of  w  are  self-explanatory  and  are  in  accord 

with  observed  phenomena  (cf.  [30]).  The  moothness  of  w  is  more  than  sufficient  to  give 

a  sense  to  (2.36)-(2.39).  The  condition  (2.42)  corresponds  to  the  pressure  condition 

(2.28).  To  see  this  one  uses  (2.43)  to  obtain 

^  1  2  12 

w  +  —  w  -  —  w 

x  2  x  2  x 

P.{f _(Vw)  -  Xw}  -  p  { - ; - %  *} 

(1  ♦  *  >  c2 


-  h  {»>.(£  -  *0  -  ^  • 


(2.45) 


On  the  dividing  streamline  w(x.|,0)  ■  y( x^ , 0)  and  the  continuity  of  the  expression  in 
(2.45)  yields  the  continuity  of  pressure  expressed  by  (2.28).  Having  w  with  the 

regularity  stated  in  part  c)  it  is  a  simple  matter  to  verify  that  (2.36)  holds  in  T+  U  T~ 

and  that  the  weak  equation  (2.39)  holds.  The  exponential  decay  described  in  part  e)  shows 
that  the  wave  is  of  essentially  finite  extent.  Of  course,  the  proof  makes  it  evident  that 

the  constant  C2  approaches  0  as  X  approaches  Xd  and  this  is  to  be  expected  in 

solitary  wave  phenomena  (cf.  [6],  [8],  [12],  [25]).  The  nonexistence  result  in  f) 
exhibits  the  delicate  nature  of  the  bifurcation  at  ( 1^,0) ,  a  behavior  quite  different 
from  that  at  a  simple  eigenvalue  wherein  the  bifurcating  branch  crosses  the  line  of 
trivial  solutions. 
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3.  Solution*  for  Smooth  Danaitiaa 

Let  pR,  n  •  1,2, ...  b*  a  sequence  of  nonincreasing  function*  in  C*[-h,1  -  h] 

which  converge  a*  n  ♦  •  to  th«  discontinuous  density  pd  given  in  (2.26).  For 

definiteness  suppose  that  p  “  p  for  -h  <  y  <  0  and  that  p  •  p.  for 

n  -  n  ▼ 

1/n  <  y  <  1  -  h.  Let  p  play  the  role  of  the  smooth  density  p  in  (2.20)  and  consider 

n  • 

solutions  (X,f),  X  «  g/c2  of  the  corresponding  problem  ( 2. 23)-( 2.25) .  Let  u  -  u(x,y) 
be  defined  by  the  relation 

c/pn(y)  u(x,y)  -  ?{y)  -  *(x,y)  .  (3.1) 

To  avoid  introducing  more  notation  we  will  refer  to  the  pairs  (X,u)  as  solutions  of 
(2.23)-(2.25).  In  an  earlier  paper  (4]  the  first  author  obtained  global  solution  branches 
for  the  flow  problem  with  a  smooth  density  and  those  results  will  be  used  here.  Hie 
description  involves  a  paraanter  yn  defined  as  follows.  Hereafter  let  I  denote  the 
interval  (-h,1  -  h)  and  define 

-  ! 

—  -  max  - - - -  .  (3.2) 

“n  vchJ(1>  /  Pn(v‘) 

For  the  smooth  density  p  the  quantity  p  plays  the  same  role  that  X.  in  (2.1)  plays 

n  n  o 

for  the  discontinuous  density.  That  is,  (p  ,0)  is  the  bifurcation  point  of  a  branch  of 

n 

solutions.  It  will  be  shown  in  Lemma  5.S  (c)  that 

lim  p„  -  X.  .  (3.3) 

n—  n  A 

Corresponding  to  e  defined  in  (2.2)  is  a  parameter  Aj,  (cf.  [4],  formula  (7.25))  which 

satisfies 

lim  . 

n+« 

For  hjj  >  0  on*  obtains  waves  of  elevation  for  the  smooth  density  pR  and  that  is  the 
case  discussed  her*.  In  Section  7.2  of  [4]  it  is  shown  that  there  is  a  maximal,  connected 
set  V _  C  *x  (Hq(8)  D  C  (I) )  of  solutions  (X,u)  containing  (p  ,0).  Furthermore,  for 


(A,u)  €  P_\{(y  # 0 ) }  th«  paraMter  X  satisfies  0  <  y  <  X  <  u  ,  where  ym  is  a 

n  n  n  n  n 

constant,  and  tho  function  u  sstisfias  u  >  0  in  S,  u(-x,y)  «  u(x,y)  in  S,  and 

3u/3x  <  0  for  x  >  0. 

For  our  purposes  the  nonlinearity  in  the  basic  equation  (1.1)  of  [4]  can  be  assumed 
smooth  and  bounded  for  a  fixed  smooth  density  Pn  (cf.  the  discussion  following  7.14  in 
[4]).  Standard  elliptic  theory  applied  to  such  an  equation  gives  local  H2  bounds  from 
local  H1  (or  L2)  bounds  and  by  embedding  theorems,  local  BSlder  estimates  (cf. 

[22)).  Hence,  if  a  sequence  of  solutions  converges  in  Hq(S),  it  will  converge  uniformly 
on  bounded  subsets  of  S,  according  to  the  Arzela-Ascoli  theorem.  Since  u  varies 
monotonically  for  x  >  0  the  argument  given  in  Lemma  4.6  of  this  paper  shows  that  u 
converges  to  zero  as  x  +  -M»,  uniformly  on  bounded  subsets  of  Hq(S).  As  a  consequence 
the  topology  of  1  x  nj(S)  is  stronger  than,  and  hence  equal  to,  that  of 
Rx  (Hq(S)  n  C°(S) )  when  restricted  to  the  set  of  solutions  Vn.  We  conclude  that  we  may 
replace  Hq(S)  O  C°(S)  by  Hq ( S)  in  the  preceding  paragraph. 

A  solution  (X,^)  of  (2.23)  gives  rise  to  a  solution  (X,w)  of  (2.36)  provided  the 
strict  inequality  ^  >  0  holds  in  S.  If  u  m  0,  then  J>(x,y)  ”  T(y),  whence 
-  1"  (y)  >  0  in  S.  Let 

Cn  _  {maximal  connected  subset  of  P„|  >0  in  S}  (3.5) 

and 

En  -  {(X,w)  |  (X,u)  e  Cn)  •  (3.6) 

This  defines  a  one  to  one  correspondence  between  Cn  and  En,  a  set  of  solutions  of 
( 2. 36)-(2. 39) .  For  nontrivial  solutions  (X,u)  we  have  3u/3x  <  0  for  x  >  0  whence 
Sw/Sx,  <  0  for  x,  >  0  by  (2.30),  (2.34),  and  (3.1).  Since  w  is  even  in  x,  and 
vanishes  at  infinity,  it  follows  that  w  >  0  in  T  if  (X,w)  e  £^{(^,0)}. 


Lemma  3.1.  The  map  from  Cn  to  En  taking  (X,u)  to  (X,w)  is  one  to  one  and 
continuous  from  R  x  Hq(S)  to  R  x  Hq(T).  For  any  rectangle  B  *  J  x  I,  j  open  (and 
possibly  unbounded)  in  R 


-rv 


■'il'wl  mk  m  l  *u*iPw.  PiuiH  nv 


' 

.•V 


,o. 

-.t' 


and 


H'(B) 


|u|  . 

H  ’(B) 


<  kJuI  . 

1  h’(B) 

<  k2|w|  . 

21  B  (B) 


(3.7) 

(3.8) 


£ 
■T .  “ 


■f-  1 


where  K1  and  K2  depend  on  n  and  K (  also  depends  on  a  lower  bound  for  ♦y.  P  being 
associated  with  u.  The  constants  JC^  and  K2  are  Independent  of  the  horlsontal 
Interval  J. 

Proof.  The  assertions  of  the  lemma  are  valid  for  any  smooth  density  and  in  the  proof 
we  write  p  for  pn>  supposing  we  have  fixed  n.  For  any  pair  (X  ,u)  e  Cn  the 
corresponding  ijiy  >  0  and  hence  the  correspondence  between  u  and  w(  via  y(x,f)  is 
well  defined  and  one-to-one.  A  further  consequence  of  $  >  0  is  the  bound 

T ( — h )  <  f  <  7(1  -  h)  in  view  of  the  boundary  conditions  (2.24).  This  bound,  together 
with  standard  elliptic  estimates  applied  to  the  seal -linear  equation  (2.23)  yield  1* 
estimates  for  derivatives  of  +  of  all  orders,  depending  only  on  p  (that  is,  on  n). 

In  particular  <  K,  where  X  depends  on  n. 

4  In  comparing  the  H1  norms  of  u  and  w  we  first  note  that  since  p  is  smooth  and 

positive  and  the  speed  c  -  (g/X)1^2  is  bounded  above  and  below  on  Cn,  it  suffices  to 
prove  inequalities  (3.7)  and  (3.8)  with  u  replaced  by 

e(x.y)  -  c/p(y)  u(x,y)  “  ?(y)  -  *(x,y)  .  (3.9) 


From  (2.20)  and  the  coordinate  relation  (2.34)  one  obtains  :{ 


It  follows  that  for  each  fixed  x  -  x1# 

/  |Vw|2dx  -  /  -  1  U»2  +  (c/p(Y(*)>  -  *  )2}dy  .  (3.12! 

I  I  *yc/p(Y<*)>  y 

For  a  wave  of  elevation  w  >  0  and  f(x,y)  <  T(y),  and  so 
p(Y(f(x,y))  >  p(Y(Y(y) ) )  -  p(y).  Thus 


y  _ 

c  /  /p(Y($(x,s) )  )ds  +  ?<-h)  -  <i(x,y) 

-h 

y  j _  -h 

>  c  /  /p ( s ) ds  +  /  c/p(s)ds  -  f(x,y) 
-h  0 


y  _ 

-  c  /  /p(s)da  -  $(x,y) 
0 


-  T(y)  -  ^(Xf y) 


(3.13; 


y  _ 

in  view  of  the  definition  (2.20) .  Since  c  f  /p(Y(f))  +  T ( — h )  -  <i(x,y)  vanishes  at 

— h 

y  -  -h,  (3.13)  can  be  combined  with  the  Poincar6  inequality  to  yield 

/  f2(x,y)dy  -  /  (T(y)  -  ♦<x,y))2 
I  I 


<  /  (/  c/p(Y(4p(x,s)))ds  +  T(-h)  -  *)2 

I  -h 


< 


,y  _ _  x2 

(/  c/p(Y(*))  +  T( -h)  -  *>)* 

-h 


-  K  /  (c/p  (Y(t)i) )  -  *  )2  • 
I  y 


(3.14 


Since 


*  <  K 

Ty 


it  follows  from  (3.12)  and  (3.14)  that  for  any  rectangle  B  -  [a, 6]  *  I 


(3.15) 


from  (3.12)  using  a  naan  value  sstimata  aa  in  (3.16).  Henca  (3.7)  holds. 

1b  sea  continuity  of  the  map  from  Cn  to  ER  consider  a  sequence  of  solutions 
(Xn,un),  m  «  1,2,...  converging  in  R  x  Hg(S)  to  (X,u)  c  Cn.  Let  ^  and  w  be  the 
respective  counterparts  of  and  u.  If  X  ■  vn  then  u  «  0  according  to  the 

description  of  Vn .  Then  since  u^  converges  to  sero  in  Hg(S),  elliptic  estimates  show 
Vu  to  be  small  uniformly  on  S  for  all  large  m.  The  y  derivatives  of  the 

A 

corresponding  $  are  bounded  below  by  a  positive  constant,  uniformly  for  a  large  and 
m 

hence  (3.7)  implies  the  convergence  of  wm  to  sero  in  Hg(S). 

If  X  <  v  ,  then  X  <  4  (X  +  u  >  <  v  for  all  large  m  and  the  function  t 
n  m  2  n  n 

corresponding  to  u  satisfies  4^  >  0  in  5.  In  fact,  there  is  a  constant  K  >  0  so 

8*  8y 

a  ^  « 

that  -j*  >  1C  on  8  and  -jy-  >  1C  on  S  for  all  sufficient  large  m.  To  see  this. 
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consider  first  the  behavior  for  large  x.  Tbe  method  of  proof  of  Theorem  2.3(b)  of  (4)  is 
easily  used  to  show  that  lujJ  <  K'e  and  |u|  <  K'e  where  B  >  K"(yn  ”  *)• 

Elliptic  estimates  show  that  and  7u  have  the  same  type  of  decay  for  large  x. 

Let  Bjj  denote  the  subset  of  S  (or  T)  where  |x|  <  N.  It  follows  from  (3.1)  and  the 
exponential  decay  that  for  some  N  >  0  and  all  m  sufficiently  large  the  y  derivatives 


of  ip  and  i|im  are  bounded  below  by  a  constant  K'  '  1  >  0.  As  noted  at  the  outset  of  the 


proof,  the  second  derivatives  of  the  function  have  a  bound  on  S,  independent  of 


3* 


n.  By  the  Arzela-Ascoli  theorem  it  can  be  assumed  that  the  functions  converge 


uniformly  to  on  B^.  Since  >0  on  B^,  it  is  clear  that  for  all  sufficiently 


large  m  there  is  a  positive  K  which  serves  on  and  on  its  compliment  as  a  lower 


bound  for  the  y  derivatives  of  i|i  and  i)i 


To  see  that  wm  converges  to  w  first  note  that,  according  to  the  lower  bound  just 
established,  (3.7)  holds  for  the  pair  (u,w)  and  for  the  pairs  (um,wm)  with  K, 
independent  of  m,  provided  m  is  sufficiently  large.  For  any  given  e  >  0  the  norm 
of  u  in  (S  -  Bjj)  will  be  less  than  e  for  N  sufficiently  large  and  the  norm  of 

u,,,  in  H^(S  -  Bn )  will  be  less  than  2e  for  m  sufficiently  large.  Then  the  norm  of 

wm  in  hVt  -  Bjj)  will  be  at  most  2K1e  and  so  it  suffices  to  show  that  wm  converges 


to  w  in  H^Bj,)  for  each  N.  Since  the  y  derivatives  of  the  have  a  positive 


lower  bound  K  for  all  large  m,  one  readily  sees  that  the  inverse  functions  y^lx,*) 
converge  uniformly  to  y(x,i|i)  on  BN  as  do  the  functions  wm  to  w  on  B^. 

Derivatives  of  the  formulas  (3.10)  and  (3.11)  provide  C2  bounds  for  the  functions  wm, 
uniformly  in  m  for  all  sufficiently  large  m  and  hence  another  use  of  the  Arzela-Ascoli 
theorem  gives  convergence  of  wm  to  w  in  C^(Bjj)  and  hence  in  H^Bjj).  q.e.d. 


Theorem  3.2.  For  each  n,  En  is  an  unbounded,  closed,  connected  subset  of  R  x  Hq(T) 
containing  (pn»0). 

Proof.  The  continuity  of  the  map  from  Cn  to  En  implies  that  En  is  connected. 


To  see  that  E_  is  closed  let  e  E  ,  m  *  1,2,..., 


converge  to  (X,w)  in 


T  X  r  «•  J  ■**  '  '  ■*  ¥  J  V1  "  !  'i* 


£ 


*>:• 


r.-v 

::v 

O' 


UVV 

* 
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Rx  hJ(T).  the  inequality  (3.8)  shows  that  the  corresponding  functions  u]|  are  uniformly 


bounded  in  H J ( S ) -  From  section  5  of  [4]  it  follows  that  the  sequence  (1  ,u  )  converges 


in  R  x  Hq(S)  to  a  solution  (X/u).  Hie  function  u  will  be  smooth  as  will  the 
corresponding  \j>,  again  from  elliptic  estimates.  In  order  to  show  that  En  is  closed  it 
suffices  to  show  that  3i|>/3y  >0  in  S. 

If  u  =  0,  then  i()  “  ¥  and  its  y  derivative  is  positive  in  S.  if  u  t  0  then 
X  <  u  and  Vu  decays  exponentially  to  zero  as  |x|  *  •  as  was  indicated  in  the 


previous  proof.  This  decay  means  that  •|i  is  positive  for  (x,y),  SXB^,  for  some  N. 


Since  p  is  constant  near  y  “  -h  and  near  y  «  1  -  h,  $  is  harmonic  near  those 
lines.  Moreover,  iji  is  nonconstant  near  those  lines  for  otherwise  it  could  not  approach 


¥  as  | x [  ♦  «•  From  the  maximum  principle  it  follows  that  is  positive  on  the 


boundary  3S  and  so  if  it  vanishes,  it  must  occur  at  a  point  (x,y)  e  S.  Suppose  this 

3» 


occurs.  Then  since  •|^  >  0  in  S,  the  gradient  of  must  vanish  at  (x,y)  with  the 


result  that 


0  <  <  K'((x  -  x)2  +  (y  -  y)2) 


(3. 19) 


1 

ca 

. 

in  a  neighborhood  of  (x,y). 

£ 
■>  i 

(x,y).  Now  apply  (3.12)  to 

containing  x.  Since  wm  is 

V  1 

integral  of  (3<i  /3y)  1  over 

IB 

S8 

Fatou’s  lemma  would  then  show 

v-1 


in  H1 
sorhood 

v-1 


The  contradiction  just  reached  also  arises  if  one  assumes  En  bounded.  For  then 


(3.8)  shows  that  C„  is  bounded  and  hence  there  is  a  sequence  (X_,u_)  converging  to 
n  n  hi 


(X,u)  e  V  \C  for  which  the  corresponding  tfi  must  vanish  somewhere.  q.e.d. 
n  n  y 


For  later  use  we  include  the  following  result* 

Lemma  3.3.  For  each  n,  the  identity  map  on  En  is  continuous  from  R  x  hJ(T)  to 
R  x  (hJ(T)  n  C  1(T)  ). 
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Proof .  The  proof  is  similar  to  the  proof  of  continuity  in  Lew  3. 1  and  we  use 


notation  from  that  proof.  In  fact  the  continuity  on  bounded  subsets  of  T  follows  from 

the  arguments  in  the  last  paragraph  of  that  proof  and  so  only  the  "tail”  behavior  need  be 

examined.  Suppose  (X _,w_)  converges  to  ( X ,w)  in  Rx  H~(T).  Then  from  (3.8)  the 
in  in  u 

corresponding  jugj  1  can  be  made  arbitrarily  small  for  large  N,  uniformly  in 


m.  Elliptic  estimates  then  show  lugj  ,  to  be 

C  (S-B„) 


11,  uniformly  in  m.  Since  the 
3* 


function  $  corresponding  to  w  satisfies  K  >  0  in  S,  one  has  t— —  >  K/2  in 

ay  ay 

S  for  the  corresponding  ^  when  m  is  sufficiently  large.  Then  the  formulas  (3.9)- 


(3.11)  show  w_|  . 

*  C  (S*Bj) 

continuous.  q.e.d. 


to  be  small  uniformly  in  m,  so  the  identity  map  is 
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4.  Sets  of  Solution*  with  Restricted  Gradient* 

In  Section  3  it  w a*  shown  that  for  «aeh  n  tha  problem 

"ET  (pn,x2)fl,7w))  "  ■Xpn(x2,w  **  T'  <4,1) 

w(xf,-h)  “  w(x1#  1  -  h)  -  0,  x1  e  R  .  (4.2) 

lim  w(x ,,x2)  -  0  (4.3) 

!*,!*• 

has  an  unbounded,  closed,  connected  set  of  solutions  En  emanating  from  the  pair 
(Mn>0).  As  stated  earlier,  the  existence  of  solutions  for  the  problem  with  discontinuous 
density  will  be  obtained  by  allowing  n  to  approach  infinity.  The  intricacies  of  the 
limiting  process  will  be  described  in  Section  5.  Here  we  describe  a  general  setting  for 


obtaining  connected  sets  of  solutions  by  letting  parasmtere  approach  limits  and  also  prove 
some  estimates  which  will  be  used  in  the  limiting  process. 

Consider  a  metric  space  X  and  a  collection  of  sets  A  -  {An>,  n  ”  1,2,...  with 
Aj j  C  X  for  each  n.  By  definition  lim  inf  A  consists  of  points  p  e  X  such  that  every 
neighborhood  of  p  contains  points  of  all  but  a  finite  number  of  the  sets  AR  while 
lim  sup  A  consists  of  those  p  t  X  such  that  every  neighborhood  of  p  contains  points 
from  infinitely  many  of  the  sets  »n.  The  following  result  from  Whyburn  ([31],  p.  15) 
provides  a  tool  for  demonstrating  the  permanence  of  connected  sets. 

—  4. 1.  let  {Ar},  n  “  1,2,...  be  a  sequence  of  connected  sets  in  a  metric  space  such 
that 

m 

(a)  U  AR  is  precompact 
n“  1 

and 

(b)  lim  inf  (A  }  f  f  . 

n 

Then  lim  sup {A  }  is  a  compact,  connected  set. 
n 

1  0. 1  - 

While  our  eventual  goal  is  an  unbounded,  connected  set  5  in  R  x  (H0(T)  n  C  (T)) 
a  first  step  will  be  the  application  of  Lenta  4.  1  to  bounded  sets  in  X  -  R  x  hJ(T).  The 
bounded  sets  to  be  considered  are  defined  as  follows.  For  S  e  (0,1)  and  R  >  0  let 
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(4.4) 


B(R)  -  {(A,w)  c  R  *  H01<T)||w|  ,  <  R  and  0  <  A  <  2Ad> 

.  H0 
and 

Qj  “  {(P1»P2>  €  «^|Pi  «  t“4  1#S  ']»  p2  «  [-1  +  8,8-1]}  >  (4.5) 

Starting  with  En  from  the  previous  section  let 

*n,5,R  “  €  En|(A,w)  e  B(R )(  Vwfx^Xj)  e  Qfi,  (x^Xj)  e  T} 

and 

F  ,  _  =  (maximal  connected  subset  of  F  .  _  containing  (y  ,0)}  .  (4.6) 

n»0#R  n,5,R  n 

Here  En  and  Fn  j  r  ue  regarded  as  subsets  of  R  x  Hg(T).  For  (A,w)  c  En,  Vw  is 

smooth  and  approaches  zero  as  |xf|  ♦  «  so  the  condition  Vw  e  is  unambiguous.  It 

will  be  shown  in  Lemma  5.5(c)  that  y  +  A  as  n  ■  so  that  for  each  pair  (6,R) 

n  d 

lim  inf{F  .  >  contains  (A_,0).  It  may  also  be  assumed  that  the  values  of  n  under 

n  r  o  t  R  d 

consideration  are  large  enough  so  that  y  <  2A .,  making  (y  ,0)  lie  in  B(R). 

n  d  n 

From  Theorem  3.2,  E  is  an  unbounded,  connected  Bet  containing  (y  ,0).  In  Section 
n  n 

5  it  will  be  shown  that  for  each  8  e  (0,1)  and  each  R  >  0  the  collection  U  {F  .  } 

n ,  o  t  R 
n 

is  precompact.  Thus  from  Lemma  4#  1 

G8,R  “  Um  suP{Fn,8,R}  (4*7) 

is  a  compact,  connected  set  containing  (A.,0).  In  Section  6  it  is  shown  that  G,  _  iB 

d  0  #  R 

still  connected  as  a  subset  of  R  x  (hJ(T)  nC0,1(f)).  Moreover,  it  is  shown  that  Gfi  R 

must  contain  a  pair  (A,w)  for  which  |w|  ^  =  R  or  for  which 

-1  "° 

max{|w  |  ,  |w  |  }  =  8  Hie  set 

x  1  »  x2  «• 

L  <T)  L  (T) 

OS 

S  *=  U  G. 

N=2  _,N 

will  then  be  seen  to  fulfill  the  requirements  of  Theorem  2.  1. 

The  restrictions  embodied  in  the  definition  of  F  ,  _  reflect  an  intimate  link 

n,6,R 

between  the  topologies  in  which  there  are  unbounded  solution  branches  and  those  in  which 
one  can  do  analysis.  The  solution  pairs  (A,w)  e  5  will  not  be  unbounded  in  R  x  C(T) 
nor  could  we  hope  to  analyze  equation  (4.1)  adequately  for  w  having  less  than  C0' 1 
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5 

« 


regularity.  The  singularity  of  f  (p  1  ,p2)  for  p2  “  -1  and  the  lack  of  uniform  convexity 

for  large  values  of  (p^.pj)  make  elliptic  estimates  elusive.  The  behavior  of 

compounds  the  difficulty.  For  solutions  in  F  .  _,  however,  one  can  obtain  the 

n,o#K 

estimates  needed  for  the  delicate  compactness  arguments  carried  out  in  the  next  section 
and  ultimately  used  in  Lemma  4.1. 

For  the  remainder  of  this  section  it  is  assumed  that  n,  4,  and  R  are  fixed  and 
estimates  are  derived  for  solutions  (X,w)  c  F  ,  We  suppress  n,  writing  p  for 

Ilf  Of  K 

p  .  The  constants  in  the  estimates  may  depend  on  one  or  more  of  the  parameters  n,  6, 
n 

and  R  and  this  will  be  indicated.  Dependence  on  p+,  p_,  and  h  is  still  suppressed. 
Recall  that  for  large  n,  the  number  X  is  bounded  above  by  2X^  for  solutions  under 
consideration.  Thus  while  the  size  of  X  enters  in  the  estimates  it  can  be  absorbed  into 
constants  depending  only  on  p+,  p_,  and  h. 


4.2.  If  ( X,w)  e  F  .  and  m  e  * 

— -  TlfOfK 


m+1  m+2  ,  m+3 

/  /  |Vwx  I  <  «S  /  /  l*r|  <  * 2  I  !  "2 

mil  *-1  I  m-2  I 


(4.9) 


where  K1  and  K2  depend  on  6.  For  e  >  0  let  I£  “  {x2  c  l| |x2|  >  e}.  Then 


m+1 

/ 

m  X 


m+2 


I  i  |Vwx  r  <  Kj  /  /  1*1 

e  2  1 


(4.10) 


where  K3  depends  on  5  and  e. 

Proof.  Let  c  e  C™(R)  be  a  function  which  has  range  in  [0,1],  has  support  in 

2 

[m  -  2,m  +  31,  and  is  1  for  x1  e  [m  -  1,m  +  2].  Multiply  equation  (4.1)  by  C  w  and 
integrate  over  T  to  obtain 


2  +  w 

1  ,,  2  _  2 

2  JJ  ***•  2 

i  t  (i  +  w  r 


,  X1  2 
1*1  -  -2  //  [pCC'  1  +  2X  //  PC  wwx  (4.11) 
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1 .■ 


after  an  integration  by  parts.  For  Vw  t  the  expression  (2  +  wx^)/( 1  +  wx^)z  is 

bounded  below  by  a  positive  constant  depending  on  6.  Thus  if  Young's  inequality  is  used 

on  the  right-hand  side  of  (4.11)  and  a  small  multiple  of  //  pc2|Vw|2  is  subtracted  from 

the  left  side,  the  second  inequality  in  (4.9)  follows. 

How  let  v  «  w_  and  f ;  j  “  - — r —  for  i,j  «  1,2.  The  result  of  applying  s— , 

*1  AJ  «Pj3Pj 

to  both  sides  of  (4.1)  is 


3 

ET7 


pfi;J(Vw) 


|j - Xp'v 

j 


(4.12) 


Let  C  denote  a  new  cutoff  function  which  is  1  on  [m,m  +  1]  and  has  support  in 

2 

[m  -  l,m  +2).  If  (4.12)  is  multiplied  by  5  v  and  an  integration  by  parts  is  used  the 
equation 


3v  3v 

3x.  3x, 
i  j 


"2  //  P«*vf  |f-  ♦  2X  //  pC2Wx 
T  j  T 


(4.13) 


is  obtained.  A  simple  calculation  shows  that  the  quadratic  form  determined  by  *i jtP-|»P2^ 
satisfies 


1  .  2  2plaia2  a2(1  *  P1)i 

ij*1'  *  “  (1  +P2>  1  1  1  +  P2  (1  4-  p2)21 

a2 

>rrp:M<1  +  7r~T2  11  ■  pi(^  ■  1)]>  <4-14> 

2  ( i  ♦  p2> 

2-2  -1 

for  all  (a.j,a2)  e  tr  and  any  n  >  0.  In  Qj,  p^  <  6  and  6  <  (1  +  p2>  <1+6 
The  choice  n_1  “  1  +  62/2  in  (4.14)  yields 


Furthermore ,  for 


f 


.  .a.a. 
ij  »  j 


+  a2} 


(PvP2*  €  °S' 


(4.15) 


(4.16) 


and  so  the  use  of  Young's  inequality,  as  before,  yields  the  first  inequality  in  (4.9). 


independently  of  n 


One  cannot  hope  to  obtain  an  inequality  Just  like  (4.9)  for  wx  , 

2 

for  the  discontinuity  in  the  limiting  density  produces  a  discontinuity  in  w  along 

*2 

x,  -  0.  However,  for  -h  <  x,  <  -  e/2,  p  is  constant  and  w  satisfies 


ST7  ri 

X 


f . (Vw)  -  0 


(4.17) 


If  w  is  extended  to  -2h  +  e/2  <  x2  <  -h  so  as  to  be  odd  with  respect  to  x2  “  -h,  it 
will  satisfy  (4.17)  in  the  extended  region.  Suppose  now  that  (  is  a  cutoff  function 
which  is  1  on  [m,m  +  1]  x  [-2h  ♦  e,-e)  and  which  has  support  on 

[m  -  1,m  +  2]  x  t-2h  +  e/2, -e/2].  Let  s(x1,x2)  m  wx2*x1'x2^'  Differentiation  of  (4.17) 
with  respect  to  x2  produces 


~  f  (Vw)  |£- 


(4.18) 


If  equation  (4.18)  is  multiplied  by  C  ■  and  the  pattern  of  the  previous  estimates 
repeated,  there  results  an  inequality  similar  to  (4.10)  but  with  the  left-hand  member 
integrated  over  -h  <  x2  4  -e.  Hie  constant  Ij  depends  on  S  and  on  e  (through  the 
derivative  |jj— )•  The  estimate  for  x2  >  e,  valid  for  all  n  such  that  p  “  Pn  is 
constant  for  x2  >  e/2  (i.e.  n  >  1/2c),  is  done  similarly.  q.e.d. 


The  next  lemma  can  be  obtained  using  standard  elliptic  theory  as  in  [26],  Lemma  3.2. 
However,  in  the  case  of  two  independent  variables,  HBlder  estimates  follow  from  the 
readily  accessible  results  of  Meyers  [18]  and  we  choose  to  use  those  here. 


Lemma  4.3.  Suppose  (X,w)  t  F  .  Then  there  is  an  a  >  0,  depending  on  6,  such 
“ — — — •  n,o,K 

that  for  any  m  e  R 


,  m+3 

i\  « ;  /  -2  • 

C  ((m,n+1]*i)  sr-2  I 


(4.19) 


lurj 

1  «  K  }  I 

*1  C®(  (m,m+1]xl)  m-2  I 


(4.20) 


\.t*'  *.’  <e";  ^*4 


and  for  n  >  1/2 e 


m+3 


l«x  I  a  *  K2  /  / 

2  C°(  [m,m+1]  *1  > 
e 


n-2  I 

where  Kg,  Kj  and  K2  depend  on  &  and  Kj  depends  on  e  as  well. 
Proof .  The  function  w  satisfies  the  equation 

h: 


(4.21) 


(4.22) 


where 


?ij<V*2>  ”  /  fij(tVw<x1.*2))dt 


(4.23) 


Since  is  convex,  f^.  satisfies  the  inequalities  (4.15)  and  (4.16)  obtained  for 

fAj.  Let  A  be  the  matrix  with  components  p?^  (i»j  “  1,2)  let  u  «  c w  where  £ 

is  a  function  of  x^  which  has  support  in  m-1<x.,<m+2  and  which  is  1  for 
m  <  x1  <  m  +  1.  Let  0  be  a  domain  with  a  boundary  30  of  class  c"  satisfying 
[m  -  1,«  ♦  2)  x  1  C  a  C  (m  -  2,m  +  3)  x  1  . 

From  (4.22)  it  follows  that 

div(A  grad  u)  -  div  f  +  g  in  0  (4.24) 

u  »  0  on  30  (4.25) 

where 

f  -  wX  grad  t.  *  (0,Xcpw) 

and 

g  *  -Acpw  +  grad  5  •  (A  grad  w)  . 

2 

From  Lemma  4.2  and  the  Sobolev  embr  ding  theorem  it  follows  that 


1*1  „  +  kl 


LP(fl) 


2  «  Kp^l  2 

L  (0)  p  L  ([m-2,m+3]xl) 


(4.26) 


for  any  p  e  (2,*).  Since  A  satisfies  a  uniform  ellipticity  condition  (cf.  (4.15), 
(4.16)),  Theorem  1  of  Meyers  [18]  yields 


t 


-30- 


| grad  uj 


(4.2?) 


<  K|w| 

LP(Q)  l/([m-2,m+3)xl) 

for  some  p  -  p(2)  >  2  and  a  K  depending  on  2  and  p(  2) .  Now  from  imbedding  theory 

(cf.  [13],  [22])  the  inequality  (4.19)  follows. 

Inequality  (4.20)  follows  in  a  similar  manner  by  ueing  equation  (4.12).  One  usee  A 

with  entries  pf^fVw)  in  place  of  A,  lets  f  »  vA  grad  C  *  (0,X(pv),  and  lets 

g  ”  -X;pv  +  grad  i  •  A  grad  v.  The  Sobolev  embedding  theorem  together  with  lemma  4.2 
2 

(with  a  cutoff  function  equal  to  1  on  Q)  yield 

|v|  <  rfw|  ,  (4.28) 

LF(fl)  L  ( [m-2,m+3]xi) 

and  Meyer's  result  gives 

| grad  v|  <  ic|w|  ,  .  (4.29) 

Lp(0)  L  [(m-2,«rt-3]*I) 

This  last  inequality  implies  (4.20)  and  will  also  be  used  in  what  follows.  The  proof  of 
(4.21)  is  done  similarly  but  with  a  cutoff  function  (  which  vanishes  for  |x2|  <  e/2. 

To  complete  the  proof  we  merely  note  that  while  the  values  of  a  occurring  in  the  three 
estimates  may  differ,  we  may  let  a  stand  for  the  smallest  one.  Likewise  we  let 
p  “  p(5)  stand  for  the  smallest  one  occurring.  q.e.d. 

The  estimates  obtained  so  far  can  be  combined  with  equation  (4.1)  to  produce  improved 

estimates  for  . 

x2 

Theorem  4.4.  Suppose  (X,w)  e  F_  ,  _  and  m  e  *•  Then 

n  f  0  t  a' 

m  «,  m+4 

(a)  (/  /  |grad{p(f  (Vw)  -  Xw)}|pr/p  <  K  /  /  w2  (4.30) 

m-1  I  m-3  I 

where  p  «  p(2)  from  Lesma  4.3  and  depend  on  2. 
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8 


|p(f,<7w)  -  Xw)T  <  K.  /  /  w2 

C  ([■,«♦  1}*I)  A  0-3  I 


(4.31) 


where  a  (as  in  Leans  4.3)  and  K2  depend  on  6. 


lW*2lc(t»,a+1)xl)  <  *3  ^_3  f 


(4.32) 


where  K-,  depends  on  S. 


(d)  Xf  p '  =0  on  a  convex  subset  Q  of  [a,a  +  1)  *  I,  then 


,  a+4 

l-x  I2  <  *4  /  /  ^ 

2ca<aj  -3  1 


(4.33) 


where  a  is  the  exponent  froa  part  b)  and  depends  on  4. 

Proof.  For  part  (a)  let  Ffx^.Xj)  ”  p(x2){fj(7w)  -  Xw>.  Recalling  the  notation 


v  -  w_  one  has 
X1 


and  from  equation  (4.1) 


+  xpV 

2  1  2 


p(f,l(Vw)v  ♦  f12(Vw)v  ♦  Xw  } 
1  2  2 


Since  | w  |p  <  «~,p~2,|w  |2  it  follows  froa  inequalities  (4.9),  (4.28),  and  (4.29)  that 

2  *2 


(4.30)  holds  with  p  and  K.  depending  on  4. 


Part  (b)  follows  froa  the  eabedding  of  W1,p  in  Ca  for  p  >  2 


For  parts  (c)  and  (d)  let  N  denote  the  L‘  norm  of  w  on  the  set 


[a  -  3, a  +  4]  x  i.  since  p  is  bounded  below  by  p+,  it  follows  froa  part  (b)  that 


|f2<Vw)  -  Xw|c( [a,»+1J*l)  *  ra 


for  a  K  depending  on  4.  Since  w  is  bounded  by  KgN  from  Lenaa  4.3,  it  follows  that 


|f2(Vw(x1fx2)|  <  K'N 


■* .  -  -  *' .  *  ,  ■  -  ■'«  *V'V 

'  ■■  ■ivwj.-vy* iv./.-.i-v a- 


»  ,  **  \\  '.S  . 


v/ 


on  [m,n  +  1]  *  I  with  1C*  depending  on  o.  Prom  the  same  le 


ww  is  likewise 
X1 


bounded.  Now 


l«2<Pl'P2)l 


2p  +  p2  -  P2 

*P2  p2  P1 


2d  +  P2> 


-1  3  ' 

2(1  +  6  )  26 


and  so  th«  estimate 


C(  [m,m+1]  *1) 


follow*. 


For  part  (d)  l.t  pt  «  wxi*x1'x2^  ,nd  “  wx  (*i'x2*  *or  *■  “  1»2  where  (x.,,x2) 


and  (Xj,x2>  are  in  fl.  From  (4.19)  and  part  (b)  the  Ca  norm  of  f2(Vw)  on  0  is 


bounded  by  KgN  where  Kg  depends  on  $.  Thus 


|f2<PvP2>  "  f2(PrP2)  +  f2<PrP2‘  "  *2(*VP2,I  *  V4*  (4,34) 

where  A  denotes  the  distance  between  (x1,x2)  and  (x^Xg).  From  Lemma  4.3  the  C( 

non  of  wx^  on  ft  is  bounded  by  KM.  Mow  for  (Pv?2)  e  Qj,  f2i  is  bounded  above 


and  f 22  is  bounded  below  by  positive  constants  depending  on  5 .  If  these  estimates  are 


used  in  (4.34),  the  result  is 


|P2  -  P2I  4  K4N4“  ' 


that  is,  w  satisfies  the  desired  HSlder  condition.  This  last  inequality,  combined 
*2 


with  part  (c)  yields  the  C  bound  in  (4.33).  q.e.d. 


Lemma  4.5.  Let  0  be  any  relatively  open  subset  of  T  on  which  the  density  p  is 


constant.  Then  for  any  solution  (\,w)  e  F_  ,  _,  w  is  real  analytic  on  0. 

n ,  o,  k 

Proof.  Extend  p  to  be  p_  for  x2  <  -h  and  to  be  p+  for  x2  >  1  -  h.  Extend 


w  to  -2h  <  x5  <  -h  to  be  odd  about  x9  “  -h  and  to  1  -  h  <  x,  <  2( 1  -  h)  to  be  odd 


about  x2  •  1  “  h.  Then  w  satisfies  (4.17)  in  each  region  where  p  is  constant. 


7  -  7  4  •  4  f  -  e  «  .  -  •  .  •  r*  .  .* 


V  ’  Vi*  ■*".  ’ 


Since  w  e  C'r~  in  ft,  the  Schauder  theory  ([13],  Chapter  6)  shows  that  w  e  C*'“(ft). 
Then  the  results  of  Morrey  ([20];  [21],  Chapter  6)  show  w  is  analytic.  Alternatively  a 
change  of  variables  back  to  4>(x,y)  which  is  harmonic  where  p  is  constant,  yields 
analyticity  for  i(i  and,  since  >  0,  for  y(x,$>)  and  w(xyX2)  (cf.  (2.34)).  q.e.d. 


The  estimates  up  to  this  point,  regarding  the  restriction  of  w  to  [m,m  +  1]  x  i, 
have  been  independent  of  m.  Estimates  which  reflect  the  decay  in  w  as  |  x  ^  |  *  «  will 
also  be  needed.  Before  giving  these  estimates  we  make  the  simple  observation  that  since 


w*  >  -1  for  any  solution  (X,w)  c 


Fn  5,R  and  8ince  w  “  0  for  *2  m  1  ”  h,  we  have 


0  <  w(x^,Xj)  <1  in  T  . 


(4.35) 


Recall  that  we  are  treating  the  case  e  >  0  (cf.  (2.41))  and  hence  w  >  0  in  T.  in  the 

jw 

case  e  <  0  one  would  have  -1  <  w  <  0  in  place  of  (4.35).  For  e  >  0,  r —  <0  on 

8x1 

(0,»)  x  I  and  this  is  used  in  obtaining  the  following  bound. 


Lemma  4.6.  For  (X,w)  e  F 


n,5,R' 


max  w(s,x»)  <  r(— ) 

v  .r  *  8 


Proof .  For  any  s  e  (0,«)  there  exists  an  x(s)  e  (s/2,8)  such  that 


_2  s 

T~  *  f  /  w  “  f  /  "2  <*»*,><**,  > 

2  s/2  I  *2  2  I  X2  22 


by  the  intermediate  value  theorem.  Since  w  is  decreasing  for  x1  >  0 


_2  x 


(4.36) 


(4.37) 


f-  >  /  /  w2  >  x  /  w2(x,x  )dx  . 


(4.38) 


Combining  these  last  two  inequalities  one  has 


v.v  V.v  v.v.v,'.  - 

«"  •“  ^  mm  t  *  -  *  •  V-  “  4  ’  a\*  "  •  *  e  ».  '  W  ■'>.'■  *V  *  1  ***  W|  »  * I* 


max  |w(x,x2)|2  <2  /  w(x,x2>wx  (x,x2)dx2 
*2eI  I  2 


Since  2x  >  s,  inequality  (4.36)  follows.  q.e.d. 


Now  we  combine  Lemma  4.3,  Theorem  4.4,  and  Lemma  4.6  to  obtain  the  following 
estimates,  wherein  a  is  the  Holder  exponent  from  the  earlier  results. 


Corollary  4.7.  Let  d(m)  «  min(1,R//m).  Then  for  (X,w)  e  F„rgrf 

M  o,1  4  K1d(B)  ' 

C  '  <[a,-)xl) 


(4. 


|w  |  <  K2d(m) 

1  c  ([m,»)xi> 


(4. 


|wx  |  <  K3d(m)  (4. 

2  ca<{t*,-)xi}ro) 

and 

|p(f_(Vw)  -  Aw) |  <  K4d(m)  (4. 

4  C°([m,-)xl) 

where  fl  is  any  relatively  open  subset  of  T  on  which  p  is  constant  and  K^» 
i  -  1,2,3,  depend  on  6,  but  are  independent  of  m,  n,  R  and  R.  The  constant  K 


1 

S 

\ 


depends  only  on  6  and  R. 


5.  Precompactne as  of  Approx xmat mg  Flows 

In  this  section  it  will  be  shown  that  U  j  R  is  precoapact  in  R  *  Hg(T) .  The 

n 

use  of  Lemma  4. 1  will  then  enable  us  to  obtain  certain  connected  sets  of  solutions  of 

problem  (4.1)-(4.3)  as  a  step  toward  obtaining  the  set  S  in  the  stronger  topology 

Rx  (Hq(T)  fl  C°'*(T))  given  in  Theorem  2.1.  Throughout  this  section  S  e  (0,1)  and 

R  >  0  will  be  fixed.  It  will  be  shown  below  in  Corollary  5.3  that  for  each  fixed 

n '  Pn  6  R  is  ptecoopact.  The  heart  of  the  compactness  question  concerns  a  sequence  of 

solutions  (1  ,w  )  e  F  ,  where  n  ■  1,2,3,...  .  The  estimates  of  the  previous  section 
n  n  n, o , R 

allow  one  to  extract  a  subsequence  of  these  pairs  which  converges  to  a  pair  (X,w),  the 
convergence  being  weak  convergence  in  Hg(T);  norm  convergence  in  C,  uniformly  on 

bounded  sets;  and  norm  convergence  in  C° ' 1 ,  uniformly  on  compact  subsets  of 

-  3 

T  -  {x2  «  0).  Moreover,  from  Learn  4.2  it  can  be  assumed  that  vn  -  wn  has 

v  -  "7"” —  w  as  its  weak  H1  (T)  limit.  Henceforth  it  will  be  assumed  that  a  subsequence 
3x1 

has  been  chosen  and  renumbered  so  that  (Xn,wn>  has  the  aforementioned  properties. 
Moreover,  it  can  be  assumed  that  wn  j*  0  for  all  n;  otherwise  there  would  be  a 
subsequence  with  second  components  all  zero  and  thus  convergent.  The  notation 
I*  -  {x2  e  l|±x2  >  0};  T*  -  R  x  I*  (from  Section  1);  Tm  -  [m,»)  x  i;  t*  -  [m,»)  x  i*, 
and  d(m)  from  Corollary  4.7  will  be  used. 


Theorem  5.1.  Suppose  (1  ,w  )  e  F  ,  „  for  n  “  1,2,...,  that  X  ♦  X,  and  that 

n  n  n, o  ,R  n 


Wn  ♦  w  weakly  in  Hg(T).  Then 


(a)  X  <  X^  and 


(b)  wn  ♦  w  strongly  in  H^B)  for  any  bounded  B  C  T. 

(c)  The  function  w  is  real  analytic  in  T  -  {x2  «  0},  even  in  x^ ,  and  satisfies 


w  -  i  <  W(m)  * 

n 


|w  I  <  Kd(m)  . 

1  C°(Tm> 


‘  *  *  '  *  *  «  •  '  »  *,  .  *„  \  ,  B  f  t  ,  * 


-.-s- 


! 

fa 


k 


Lemma  5.2.  Suppose  (X,w)  €  F_  ,  „  and  that  0  *  1  -  (X/u_)  >  0.  Then  for  m  >  0 

n  i  o  ,  n  n 


M  i  < 
Hn<Tm) 


where  K  depends  on  5 ,  R,  and  8 . 


Proof.  Since  |w|  <  <  R  by  assumption#  the  inectualitv  (5.6)  need  be  demonstrated 

h’<t> 

only  for  m  larger  than  some  mg.  If  equation  (4.1)  is  multiplied  by  w  and  integrated 
over  it  yields 


7  /  /  P  - \  |Vw|  2  =  -  /  p  w(m,x  )f  (Vw)dx  -  X  /  /  p'w2 

2  T  n  <1  +  w  )2  n  21  2  Tn 

m  x  m 

2 

<  -  !  Pnw(m,x2)f 1(Vw)dx2  +  L-  /  /  Pn | Vw | 2  (5.7) 

I  *n  T 

m 

from  the  characterization  (3.2)  of  y^ .  For  m  >  m0,  mg  depending  on  6  and  R, 

(2  +  w„  )/2(1  +  w  )2  >  1  -  (0/2)  and  -wf.(Vw)  <  2ww„  for  x.  >  m.  Hence,  from  (5.7) 

x2  x2  1  xi  i 

j  0P+  /  /  |Vw|2  <  p_  /  |w(m,x2) | |w  (m,x2)dx2| 

T  I  1 

m 

which,  with  Corollary  4.7  and  the  Poincar6  inequality,  provides  (5.6).  q.e.d. 


Corollary  5.3.  For  each  fixed  n'  Fn  $  r  is  PrecomP*ct  in  R  x  Hq ( T ) . 

Proof.  Let  (X,.  k  »  1,2,3,...,  be  a  sequence  in  F_  *  If  a  subsequence 

■  x  x  n,o,K 

converges  to  (y  ,0),  we  are  done.  If  not,  then  | X,  -  y  j  +  |w.|  .  >  n  >  0  for  all 

n  x  n  K  „  ■ , _ . 

H  (T) 

large  k.  Suppose  a  subsequence,  still  denoted  ( X^ ,wk )  is  such  that  X^  converges  to 
Un  as  k  +  Since  wk  lies  in  the  ball  of  radius  R  in  H ^ (T) ,  u^  lies  in  the  ball 
of  radius  K2R  in  H^S)  by  (3.8).  According  to  [4),  section  5,  the  collection 


{(X^.u^)},  k  •  1,2,...  is  precompact  in  ft  x  (L‘(S)  D  Cq(S))  (the  exponent  o 


in  1 . 4  of 


[41  is  1  as  established  in  the  discussion  following  (7.14)  of  that  paper).  It  follows 
that  a  subsequence,  still  denoted  (X^jU^),  converges  to  (wn,u)  in 

R  *  (L2(S)  O  Cq ( S ) ) .  From  local  elliptic  estimates  one  sees  that  (yn,u)  is  a  solution 


>*AV-  >**'.*■  • .  •  .*•  ■.■> •. ‘A4-  .*»  1--  *•-  .'•'•>  VA,'-  vAv  AV'  aV>Y 


m®, 


2  - 

and  that  tha  uk  converge  to  u  in  CQ(S) •  Then  Theorems  2.3  and  2.4  of  [4]  show  that 

u  “  0.  From  this  information  ons  can  conclude  that  ths  y  derivatives  of  the 

corresponding  ^  are  bounded  below  by  K  >  0,  uniformly  on  S  and  uniformly  in  m 

for  m  large.  Elliptic  estimates,  uaed  once  more,  show  that  the  L  convergence  of 

uk  to  aero  implies  convergence  to  sero  in  H 1  (8) ,  and  (3.7)  implies  convergence  of  wk 

to  sero  in  H^T).  However,  from  an  earlier  inequality  it  follows  that 

|wv|  ,  >  n/2  >  0  for  all  large  k,  a  contradiction. 

*  h’(T) 

From  the  last  paragraph  it  is  clear  that  X^/i^  <  6  <  1  for  all  k  and  thus  (5.6) 
holds  for  all  w^,  uniformly  in  k.  The  precompactness  of  {wk>,  k  ■  1,2,...,  in 
Hg(B)  for  each  bounded  set  B  C  T  follows  as  in  the  proof  of  Theorem  5. 1.  These  two 
properties  combine  to  yield  precospactness  in  *  x  hJ(T).  q.e.d. 

To  show  precompactness  of  U  Fn  ^  R  it  will  now  suffice  to  show  that  the  sequence 

n 

{(X  , w  )},  n  •  1,2,...,  selected  at  the  outset  of  this  section,  converges  in  Rx  h,!(T). 
n  n  w 

He  next  show  this  together  with  some  of  the  properties  of  the  branch  S  required  for 

Theorem  2. 1.  For  the  case  lim  A  -  A  <  X  .,  we  need  merely  quote  results  already 

n  & 

established  to  see  the  convergence  of  wR.  For  X  -  Xd,  a  delicate  analysis  must  be 
carried  out  using  the  nonlinearity  in  equation  (4.1). 

The  case  X  <  X^ 

Theorem  5.4.  Suppose  (*n#"n>  <  Fn,«,R  satisfy  Sj  ^  X  4  Xd  “d  wn  +  in 

Hg(T)  as  n  ♦  “•  Then 

(a)  wR  converges  to  w  strongly  in  Hg(T) . 

dw 

(b)  w  >  0  in  T  and  t —  <  0  for  xi  >  0  and  X2  e  1. 

®X1 

Proof,  (a)  It  will  be  shown  in  Lemma  5.5(c)  that  uR  approaches  Ad  as  n  ♦  "• 
Hence  inequality  (5.6)  holds  for  wn,  uniformly  in  n  for  all  large  n.  This  decay  in 
the  tail  together  with  Theorem  5.1  (b)  yield  strong  convergence  in  Hg(T). 


(b)  Pro*  (5.S)  it  follows  that 


for  all  nonnegative  ?  c  C^(T) ,  where  f^  is  defined  by  (4.23).  Hence  from  the  strong 
aaximum  principle  ((13],  theorem  8.19),  either  w  =  0  or  w  >  0  in  T.  the  function 
v  ~~  is  in  Hg(T)  by  Tamms  4.2,  and  is  the  weak  H1  limit  of  vn  -  wR.  The 
equation  (4.12)  is  satisfied  by  (\n,wn,vn)  and  since  vn  <  0  for  x1  >  0  one  obtains 

//  pdfijtVw)  4  0 


for  all  nonnegative  »  e  CQ(  {<),•)  x  I).  From  the  strong  maximum  principle  v  s  0  or 
v  <  0  on  (0,<*)  x  I,  corresponding  to  the  case  w  s  0  or  w  >  0,  respectively.  If 
w  =  0,  then  the  functions  *n,  which  are  even  in  x  ^  and  decreasing  in  xj  on 
(0,«)  x  1,  aust  converge  uniformly  to  zero  on  T.  The  gradient  estimates  (4.20)  and 
(4.32)  show  that  VwR  converges  uniformly  to  zero  on  T  and  thus  wn  converges  to  zero 
in  C0'1^).  The  analogue  of  (5.7)  taken  over  all  of  T  yields 


2  +  i 

_  - 2  |vw  !2  <  -2.  //  p  |Vw  |2 

T  n  2(1  ♦  (w  )  )2  "  vn  n  n 


//Pm 


(5.8) 


Since  -jj-  converges  to  —  <  1  and  |  Vw^J  ♦  0  uniformly  on  T  as  n  ♦  •,  (5.8)  leads 
to  a  contradiction  for  large  n.  Thus  w  >  0  on  T  and  -jj—  <  0  for  Xj  >  0.  q.e.d. 


The  case  X  ■  A . 

1  d 

As  stated  earlier,  it  ia  assumed  that  we  are  working  with  a  sequence 
(An,wn)  €  Fn ^ R  for  which  Xn  ♦  X  ■  Xfl  and  wn  converges  to  w  in  the  weak  topology 
of  hJ(T)  as  n  ♦  •.  In  addition,  the  conclusions  of  Theorem  5.  1  hold.  The  plan  here  is 
to  show  in  Theorem  5.6  that  the  weak  limit  is  w  =  0  and  then  to  show  in  Theorem  5.7  that 
wn  ♦  0  strongly  in  Hg(T)  so  that  the  limiting  point  is  the  bifurcation  point  (1^,0). 
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« 


To  begin  recall  that  un  is  defined  by 


•n2 


-  / 


'»  ucHq<z)  /  Pn«u’> 


2  * 


(5.9; 


2 

Let  9  *  0  (x- )  be  the  associated  positive  maximizer  normalized  by  /  p  9  “  1* 

n  n  2  j  n  n 


The 


function  9  satisfies  the  Euler  equation 
n 


(p  0  * )  *  -  u  p*9  . 

lHn  n'  “nHn  n 


(5.  10 


An  associated  quantity  aR  is  defined  by 


a  (X)  - 
n  t 

uchJ(I) 


-X  /  p'u  -  /  p  (U* ) 
I  n  I  n 


.,2  1 


{ 


(5.  11 


for  each  X  e  *.  Let  y„  “  Yn(*2*  be  the  positive  maximizer  in  (5. 11)  when  X  ” 


V 


aqain  suppose  that  f  P  Y  “  1«  The  function  Y  satisfies  the  Buler  equation 
•  "  j  n  n  n 


Ipj' ) 1  -  X  ft'v  -  <*  (X  >p  Y  . 
nn  n  n  n  nnnn 


(5.  12 


Of  course,  0  and  y  both  vanish  at  the  endpoints  of  I  "  (-h,1  -  h). 
n  n 


Lemma  5.5.  Let  (Xn,w„)  be  in  *nf8fR*  «»•»» 


(a)  o  ( X  )  <  0 
n  n 


(b)  Assume  Xn  ♦  Xfl  as  n  +  •.  Then  the  functions  {9n)  and  {y„} «  n  “  1,2, 

are  bounded  in  c’u)  uniformly  in  n  and  for  each  e  >  0  converge  in 
.1,5 


C(f)  n  c'(ie),  I(  ■  IH  {|*2I  >  the  function 


a(x2  +  h) 


x2  <  0 


ah 


(5. 13 


<x,  -  1  +  h) 


1  '*2 


x2  >  0 


where  a2  “  3/h2(p_h  +  p+(1  "  h)).  Let 


Q<u)  - 


”  i  P(*U  <P_  ”  P+»«2<°) 


W*')2  /Pd‘u*>2 


(5.  14 


Then  yA  is  the  positive  maximiser  of  Q(u)  for  u  e  Hq(I)  and 
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whore  X^  is  defined  by  (2.40) . 

(c)  lim  u  “  X  . 

n  a 

n+« 

Proof,  (a)  From  the  characterization  (5.11)  is  a  strictly  increasing  function 

of  X  and  from  (5.9),  a  (y  )  *  0.  As  noted  in  the  proof  of  Corollary  5.3,  for 

n  n 

(X  ,w  )  e  F„  *  „\{(u_»0)},  X  <  u  »  Mi  SO  a  (X  )  <  0. 
n  n  n  ,0  ,k  n  n  n  n  n 

(b)-(c)  Since  evaluation  at  x2  “  0  is  a  compact  mapping  for  functions  in  hJ(I), 
the  maximum  of  Q(u)  is  attained,  and  the  form  of  Q  ensures  that  this  maximizer  must  be 

linear  on  I*  and  I~.  K  calculation  shows  that  fd  is  the  positive  maximizer  which  has 

2  “1 
the  normalization  /  P^T ^  “  1  and  that  Q (Yj)  “  ^d  • 

Suppose  the  maximizer  6  for  (5.9)  is  renormalized  to  have  max  0  -  1  on  I  and 

n  n 

that  its  maximum  occurs  at  x  •  Since  0‘(x,)  >0,  an  integration  of  (5.10)  from  x 

a  n  2  « 

to  x2  yields 


pn0A(x2>  “  /  pA0n 


and  so  |0^(x2)|  «  (p_  -  p+  )/p+.  This  bound  easily  implies  that  the  eigenfunction 

normalized  by  f  p  02  »  1  must  also  satisfy  |0'|  <  K.  where  K.  depends  only  on  p. 

j  n  n  n  •  *  x 

and  h.  On  I£,  that  is,  where  |x2|  >  e,  the  derivative  -  0  for  n  >  1/e  and 

hence  6"  -  0.  It  follows  from  the  Arzela-Aacoli  theorem  that  a  subsequence  of  {0  } 
n  n 

converges  as  asserted  to  a  function  0^  which  must  be  linear  on  I+  and  I  •  Since 
/  p  02  -  1,  a  calculation  shows  that  0^  must  be  yd  and  the  uniqueness  of  the  limit 
shows  that  the  whole  sequence  converges.  The  maximum  value  in  (5.9)  converges  to 


(p_  -  p+)yd(0) 

f  °d(^)2 


and  thus  p_  ♦  X.  as  n  ♦  “. 
n  a 


All  values  of  X  under  consideration  for  a  fixed  n  lie  in  the  interval  (0,un) 

and  since  a  is  increasing,  a  (0)  <  a  (X )  <  a  ( p_ )  “  0.  Moreover,  the  quantity  a  (0) 
n  n  n  n  n  n 

is  readily  seen  to  be  bounded  below  by  a  negative  constant  K2  depending  on  p^  and  h, 
but  not  on  n.  By  integrating  the  equation  (5.12)  for  fn  •»  was  done  with  (5.10)  one 
finds  that  |y'|  <  K,  on  I  and  |t“|  <  K.la  <X_)|  on  I  for  n  >  1/e.  To  see  that 
Yn,  like  0n,  converges  to  Yd<  it  suffices  to  show  that  | aR ( Xn )  |  +0  as  n  +  «. 

Since  X  X.,  u  ♦  X.,  and  a_(y_)  ”  0,  the  desired  convergence  will  follow  if  it  is 

n  a  n  a  n  II 

shown  that  a  has  a  modulus  of  continuity  independent  of  n.  To  see  this  let  u_  be  a 
n  * 

maximizer  associated  with  <*n(r)  in  (5.11)  and  suppose  /  Pn“r  “  Then  for  r  >  s, 
an(r)  -  -r  f  p;«2  -  J  Pn(u;)2 

-  -(r  -  •)  /  p;u2  -  .  /  P;«J  -  f  pn(u;>2 

<  (r  -  s)K_(p  -  p  )  +  a  (s)  (5.15) 

d  +  n 

since  ur  is  an  admissible  function  in  the  definition  of  an(s).  Since  <*n  is  an 
increasing  function,  it  follows  from  (5.15)  that 

lan(r)  -  ®n<a>l  *  *6 I r  ”  *1 

with  Kg  independent  of  n.  q.e.d. 

Theorem  5.6.  let  S  e  (0,1)  and  R  >  0  be  fixed  and  suppose  «  Fn,S  ,R  tor 

n  ■  1,2,...  .  If  X  +  X .  and  w_  converges  weakly  to  w  e  H«(T)  as  n  ♦  •,  then 
n  <2  ** 

(a)  w  =  0  in  T. 

(b)  | wn |  and  |Vwn|  converge  to  zero  uniformly  in  T  as  n  +  •• 

Proof,  (a)  The  argument  revolves  about  the  behavior  of  wn  for  large  Xj  and  that 
behavior,  as  we  shall  see,  is  reflected  in  the  projection  of  wn  along  the  eigenfunction 
Yn  of  (5.12).  Let 


V*1'x2>  “  WW  +  Vx1'x2> 


(5. 16) 


With  the  um  of  (5.16)  tha  second  integral  In  (5.24)  can  ba  written 


4  /  T'P  (Fl'MM  )2  -  |  r2  I  p  (f* ); 

2  i  n  n  n  n  n  Xj  2  n  v  n  n 


*  *  f  l  wvL 

n  2  x  2 


(5.25) 


2  2 

where  hara  and  in  what  follows  tha  convention  (Rh)Jt^  is  used  for  { ( )  .  Mow  use 
tha  definition  of  y .  in  (5.13)  to  cosjpute 

<3 


•  “  llM  i  {  pn(Yn)3  *  T  /  pd(Yd)3  "  T  f  ~  T~  ~~^T37?  ft  "  *  (5*26) 

n-»"  1  1  [p_h  +  p+(1  -  h)]  h  (1  -  h) 

2  2 

Sinoe  the  expression  (p_/h  )  -  (p+/( 1  -  h)  )  -  e  (cf.  2.2)  is  positive  by  assuaptlon, 
a  >  0.  Only  the  positivity  of  a  is  iaportant  in  what  follows  so  to  simplify  notation 
we  ondt  the  tilde.  We  conclude  that  for  any  c  >  0  and  all  n  >  N(e) 


The  use  of  this  last  inequality  in  (5.25)  together  with  Young's  inequality  gives 

ll  !  WVi,  -  "S'  <  2**J  *  'l  !  WL  >5'”> 


where  depends  on  e ,  but  is  independent  of  n  and  x^.  From  the  form  of  A  given 

in  (5.23)  and  fro#  the  decay  of  wB  and  Fwft  as  x1  ♦  given  by  (4.39),  it  follows 
that 

|A(7wn)|  <  Kj((F^)2  +  F2  +  | VR„| 2 )  (5.28) 

for  Xf  larger  than  loae  value  X.  The  use  of  the  last  two  estieates  in  (5.24)  yields 


|g'  +  a  (X  )P  +  .f2|  <  2er‘  +  x(e)  /  p  {(F’)2  f  F2  +  |vr|2}  (5.29) 

'  n  n  n  n  n  n  j  n  n  n 


for  n  >  N(e)  and  x1  >  X  where  *(e)  is  independent  of  n  and 

Note  that  the  term  Gn  approximates  F^  ( x  j )  for  wx^  small  and  thus  (5.29)  is  akin 

to  a  differential  inequality  for  Fn>  The  aim  is  to  show  that  in  the  limit  as  n  +  •  and 

o  (X  )  +  0,  (5.29)  is  inconsistent  with  having  a  positive  function  w(x.,x,)  which 

n  n  i  * 

decreases  in  x  1  for  x  1  >  0.  To  arrive  at  this  contradiction  we  first  estimate  the 
integral  in  ( 5. 29 ) . 

If  the  principal  equation  (4.  1)  is  multiplied  by  wR  and  the  two  sides  are 

integrated  over  I  there  results 

/  p  (w  )  f_(7w  )  -  -1  /  p’w2  +  /  p  w  x— —  f  ,(Vw  ) 

Kn  n  Xj  2  n  n  ^  n  n  ^  n  n  dx^  1  n 


-  -X  /  p’w2  f  /  p  w  f  ,(Vw  )  -  /  p  (w  )  f  ,(Vw  ) 
n  '  n  n  dx  1  •  Kn  n  1  n  ^  Kn  n  x  1  1  n 


(5.30) 


°n<*1>  “  /  PnVl(7V 

x  4.  4-  n 


(5.31) 


Then  since  f2tPl'P2*  "  P2  "  (2P2  +  +  Pi)/2M  +  P2'  '  ®<Iuation  (5.30)  yields 


n  x 

/  p  (w  )2  +  /  p  — — j —  »  -X  /  p’w2  f  U  ( x . ) 

L  n  n  x„  i,  n  1  ♦  (w  )  ninn  dx  nl 

x  •  x  n  x  ^  •  ( 

+  i  /  P„  - ~—~~2  {*<"n>x  +  3(wn)2  f  (w  )2  } 

2  I  "(If  (w  )  )2  "  x2  "  *2  "  X1 

n  Xj 


(5.32) 


Now  from  the  representation  wR  -  Fnyft  f  R„,  the  equation  (5.12)  for  yn,  and  the 


orthogonality  (5.19) 


/  p  (w  )2  -  F2  /  p  (Y  ’  )^  -  2X  F  J  p’y  R  f  /  p  (R  )“ 

£  n  n  Xj  n  ^  n  n  n  n  ^  n’n  n  ^  n  n  x 


A  second  use  of  the  representation  for  w  yields 


f  p  (w  )2  f  X  /  p’w2  -  -a  (X  )F2  f  X  J  p’R2  f  /  p(R  )2 
|  n  n  n  '  n  n  nnn  n  £  n  n  ^  n* 


(5.33) 


(5.34) 


v  ■  i  o  ,i>.\  -.v.v 


'  '►  ■'  .v  .  ,*•  • 


*  ’•  «  •  .  -  -  .  ■ 

.  *  *  •  »  ■  . 
-  1  d  *  ,  •<  ,  *  .  « 


■"  w*.'.  ■  7  V  ipw  *  «p_  j\ .  ip.  i.-  i 


Now  define  t  by 
n 


i.  - ;  * 

/  “Pjr  “° 
j  n  n 


(5.35) 


Since  an  extra  condition  on  a  has  bean  added  to  the  characterisation  of  u  in  (5.9)  it 

n 

is  natural  that  t  >  u  and  this  is  shown  in  ((4),  Lena  3.1).  An  argument  similar  to 
n  n 

that  given  in  Leama  5.5  shows  that  as  n  ♦  •  the  numbers  tr  approaches  a  limit  T„  > 

X,  so  that  X  /t  +  X ,/x_  <  1.  Thus  1  -  X  /r  >  6  >  0  for  all  large  n.  Since 
n  n  a  •  n  n 

/  Vx1'x2,Vx2>Vx2,dx2  "  0 


for  each  x .. 


/  K*l  <  -r  { wL  «"-•>/  *.<*.»; 

X  n  x  «  a 


The  use  of  this  last  inequality  in  (5.34)  gives 

.2  .  .  r  ,.2 


/  *  »„  /  k<  >  -w;  *  •  /  w; 


(5.36) 


for  all  large  n.  The  inequality  (5.36)  will  be  used  in  (5.32).  First,  however,  note 
that  because  of  the  decay  given  by  Corollary  4.7,  for  any  e  >  0  there  is  an  X(e)» 
independent  of  n,  such  that 

,2 

i  o 

(5.37) 


<*  >* 
n  x 


[  Pn  1  T  (717  *  (1  _e>  l  Pn<W">xl 


provided  x,  >  X(e).  This  same  decay  coupled  with  the  decomposition  (5.16)  of  wn  and 
Young's  inequality  gives 
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‘0  jl*  -  jk’ *  j,*  * 

jd  i\lM 


t"1 

t-i 

cO 

3 

•  * 

'yi 


'a 

3 


4  w 

N  ' 

‘.H 


3 


a 


7  /  P 


2i  “  (!♦(.»  )2 

n 


n  %2 

,3.„_  r  _  ,2 


{2(w  )  +  3(w  )  «•  (w  )  } 


n  x_  n  x. 


<  n  r  +  K,e  /  p_<R  )"  ♦  c  /  p  (w  )“ 
i  n  *  ^  a  n  x  ^  j  n  n  x  ^ 


(5.38) 


for  x,  >  X(e)  (taking  X(e)  larger  if  necessary).  The  last  three  inequalities,  with 
e  chosen  so  that  K^e  <  8/2,  combined  with  equation  (5.32),  produce  the  inequality 

T  l  Pn(Rn>x,  -  an(Xn,F^  +  (1  “  2e  >  [  pn(Vx.  4  °A<*1>  +  K1*n  <5*39> 


valid  for  n  >  N(e)  and  Xj  >  X(e).  To  further  decaapose  terms  using  FR  and  Rjj  note 
that 

/  P  (w  )2  -  (F* )2  /  p  Y2  ♦  2F'  /  p  y  (R  >  ♦  /  p  (R  ) 2 

'  n  n  x.  n  '  n  n  n  j  n  n  n  x^  '  n  n  x1 


(f;)2  +  /  P(R  )2 

XI  j  X\  X  ^ 


(5.40) 


using  (5.18)  and  (5.19).  It  can  be  asaumd  that  e  is  chosen  so  that  8/2  <  1  -  2 e. 


Then  since  a  (X  )  <  0  it  follows  frca  (5.39)  and  (5.40)  that 
n  n 

|2  „  2 


<FA>2  ♦  /  pnl ™nr  <  f  °a  *  • 


(5.41) 


For  these  same  parameter  ranges,  the  use  of  (5.41)  in  (5.29)  yields 

|g*  ♦  a  (X  )F  eF2|  <  2eF2  +  K0'  +  XF3  (5.42) 

'nnnnn  nnn 

where  e  is  positive  and  K  depends  on  e ,  but  not  on  n  or  xj.  Since  Fn ( x 1 )  +  0 
as  xi  •,  uniformly  in  n,  we  may  assume  that  X(e)  is  chosen  so  that  KFn(xi)  <  e 
for  XJ  >  X(e),  leading  to 

|G«  +  o  (X  )F  +  eF2|  <  3eF2  +  KO'  .  (5.43) 

'nnnnn'  n  n 

Now  one  can  integrate  the  inequality  (5.43)  (without  the  absolute  value  taken  on  the  left 

a  A 

side)  over  the  interval  [x^x.,]  where  X(e)  <  x(  <  x^  <  •.  Next  let  n  +  <»  using  the 
results  of  Section  4  and  the  fact  that  a_(X_)  ♦  0.  Finally  let  x,  ♦  •  to  obtain 


-  /  pdYa  \  *  V"  <  -<•  -  3*>  /  -  %  J  Pdwf  ,(7*) 

1  x2  *1  1 


(5.44) 


Hera  w  is  tha  weak  limit  in  H«(T)  of  ths  aaquanca  {w  }  and  F  ia  dafinad  in  analogy 


with  Fn  in  (5.7).  Inequality  (5.44)  can  ba  writtan  as 


*»(x  >x,)  , 

lx  d  2  x, 


(5.45) 


Now  from  tha  naan  valua  theorem 


w(*l<*» 


<  i|w  (X  1# •  >  |  „ 


(5.46) 


and  front  Theorem  5.1(c)  tha  right-hand  aida  of  (5.46)  approaches  saro  aa  x1  ♦  ».  It 


followa  front  tha  last  two  inaqualitias  that  for  all  larga  x1 


2  I  pdv«  Yd  >  <•  -  3e)  /  F2(a)ds 
I  I  x, 


(5.47) 


Aa  in  tha  proof  of  Theorem  5.4  ona  concludas  front  tha  strong  naximust  principle  that 


w  <  0  in  (0,-)  x  I  or  that  w  =  0  on  T.  Tha  inaquality  (5.47)  insures  that  only 
X1 


tha  lattar  can  occur  completing  tha  proof  of  part  (a)  of  Theorem  5. 1. 


For  part  (b)  recall  that  wn  was  taken  to  converge  to  w  uniformly  on  compact 


sets.  Since  w_  has  its  maximum  on  the  line  x.  »  0,  the  convergence  to  w  =  0  is 
n  1 


uniform  on  T. 


4.3  and  Theorem  4.4  show  that  |Vw|  ♦  0  uniformly  in  T  as 

n 


n  +  q.a.d. 


Theorem  5.7.  Let  S  t  (0,1)  and  R  >  0  ba  fixed  and  suppose  U„*wn)  «  Fn,S,R  Cor 
n  -  1,2,...  .  Assume  XR  ♦  Xd  and  w„  converges  weakly  to  taro  in  hJ(T)  as  n  *  - 


Than  wn  convargas  strongly  to  taro  in  Hq(T). 


Ivi  v"  v'v  1.1 


Proof .  In  the  previous  proof  it  was  shown  that  the  inequality  (5.43)  holds  for 


n  >  N(e)  and  >  X(e),  the  latter  restriction  being  made  solely  to  insure  that  wn 


,0,1, 


sufficiently  small  in  C  *  I).  From  the  previous  theorem  we  know  that  wn 

0,1 


approaches  zero  in  C  *  (T)  and  thus  inequality  (5.43)  can  be  shown  to  hold  for  all 
x1  >  0.  This  will  enable  us  to  relate  the  Hq  norm  of  wn  to  the  quantity  Fn(0)  in 
such  a  way  as  to  show  that  wn  approaches  zero  strongly  in  hJ(T)  as  Fn(°)  ♦  0. 


From  (5.43)  one  obtains  the  two  differential  inequalities 


{G  -  KU  }  <  -a  (X  )F  -  (e  -  3e)F2 
dx1n  n  nnn  n 


(5.47a) 


and 


{G  +  KU  }  >  -a  ( X  ) p  -  (e  +  3e)F* 
ax  ^  n  n  nnn  n 


(5.47b) 


valid  for  n  >  N(e)  and  for  all  Xj  >  0.  As  noted  earlier,  we  have 


(w  ) 
n  x . 


Kw 


^  1 ♦  (w r  *-r> 

I  n  x,  n 


(5.48) 


Since 


Wn<x  1'*  > 


V*> 


<  K  Vw 


(5.49) 


L  (I) 


and  Vw 


'L  (I) 

♦  0  as  n  ♦  ■*,  N(c)  may  be  increased,  if  necessary,  so  that  for  all 


L  (T) 


n  >  N(e  ) 


(1  +  e)F'  -  (1  +  e)  /  p  Y  (w  )  <  G„  ±  KU 

n  I  n  n  n  X1  n  n 


<  (  1  -  c )  /  p  Y  (w  ) 
i  Kn'n  n  > 


(1  -  OF' 


(5.50) 


Recall  that  FI  is  negative. 


The  expression  o  (X  )  +  (e  +  3e)F  fx.)  is  monotone  decreasing  on  0  <  x,  <  <»  and 
n  n  n  i  » 

approaches  an<Xn)  <  0  as  x1  ♦  •.  Let  [Zn,«)  be  the  subset  of  [0,«)  on  which  the 
expression  is  nonpositive.  Suppose  first  that  Zn  -  0.  If  the  inequality  (5.47b)  is 
multiplied  by  the  positive  quantity  -(Gn  +  KUn)»  t*le  result  Integrated  from  x^  to  », 
and  use  is  made  of  (5.50),  then 

<1  +  e>2(FM2  >  -<1  -  e)(oU  )F2  +  4  (e  +  3e)F2)  .  (5.51) 

n  n  n  n  J  n 

Since  F ■  ( 0 )  *  0,  it  follows  from  (5.51)  that 
n 

F2(0)(a  (X  )  4(tt  3e)F  (0))  >  0  .  (5.52) 

n  n  n  j  n 

However,  the  assumption  that  Z  »  0  or  that  a  (X  )  +  (e  +  3e )F  (0)  <  0  is  in 

**  n  n  n 

contradiction  to  (5.52),  given  that  a  (X  )  <  0. 

n  n 

We  now  know  that  Zn  >  0  and  that 

o  (X  )  +  4  (e  +  3e)F  <Z  )  <  a  (X  )  +  (e  +  3e)F  (Z  )  <  0  .  (5.53) 

n  n  3  nnnn  nn 

Using  this  inequality  in  (5.51)  one  finds 

(1  +  e)2(FMx1))2  >  |  (1  -  e )  (e  +  3e)F2(Xl)<Fn<Zn)  -  Fn(x,))  (5.54) 

on  [Z  ,«•).  This  last  inequality  yields 
n 


<  K1Fn/2<V 

<  K,F*/2(0)  .  (5.55) 

On  the  Interval  (0,Zn) 

-a  (X  )  <  (e  +  3e)F  (x,) 
n  n  n  * 

and  the  use  of  this  in  (5.47a)  yields 

—  (G  -  ru  )  <  ((e  +  3e)F2  -  (e  -  3e)F2) 
dx  j  n  n  n  n 

-  6EF2  .  (5.56) 

n 

We  now  carry  out  a  process  similar  to  the  one  above  of  multiplying  (5.56)  by  ”(Gn  ”  KUn># 
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integrating  from  0  to  Xj,  and  using  (5.50)  to  obtain 


(G  -  KU  )2  | x  _  P3  ix, 

-2— _JL_  ^ecdte)-!  \ 

2  l0  3  >0 


x^  e  [0,Zn]  • 


(5.57) 


Now  (5.50)  can  be  used  to  show  that  Gjj  -  KUn  >0  at  Xj  “  0  and  to  bound 
(G^xj)  -  KOn(x1))2  above  by  (1  +  e)2(P^(x1))2.  Hence  from  (5.57)  one  has 
+  e)2(F'l 


(1  +  e)‘(F'(xJ)2  >  4e<  1  +  e) (F3(0)  -  F*(x.>)  . 
n  i  n  n  1 

valid  for  xf  e  [0,ZnJ.  He  use  this  estimate  as  follows: 


(5.58) 


Z«  Fn(0)  2 

tv2,  \A  t  t2dt 

/  Fn  X1  ^*1  J  -F ' 

0  n  F  (Z  )  n 

n  n 


<  *•  / 


F  ( 0 )  - 

n  t2dt 


VV  ^F3(0,Tt3 


and  combining  this  with  (5.55)  yields 


<  K2fV2(0) 


|  F^x^dx,  <  K3F3/2(0) 


(5.59) 


Now  we  use  inequality  (5.41)  which,  like  the  inequality  (5.43)  used  at  the  outset  of 
the  proof,  will  hold  on  0  <  x1  <  -  provided  n  is  sufficiently  large  thereby 
guaranteeing  that  wn  is  sufficiently  small  in  C^' 1  (T) .  From  (5.41)  it  follows  that 


/ 

0 


<  K, 


Fn**1 


(5.60) 


since  /  U^dx1  -  0.  the  Poincar<  inequality  yields 
0 


[  l  pn*2n  <*4  [  l  pnt7Rn|2 


(5e61) 


% 


How,  sine*  wn  1*  a  van  in  Xj  and  haa  the  fora  wn  ■  rntn  ♦  *j, 

followa  froa  the  laet  three  Inequalities  that 

|wj2,  <  K(F2/2(0)  +  f!*/2<0)>  . 

i  n'Hl(T)  n  n 


(cf.  (5. 16)-(5.  19))  it 


(5.62) 


By  Theorea  5.6  Fn(0)  ♦  0  as  n  ♦  »  and  hence  wn  converge*  strongly  to  sero  in 

hJ(T).  q.e.d. 


the  result  just  completed,  together  with  Corollary  5.3  and  Theorea  5.4,  yields 

Theorem  5.8.  For  each  6  e  (0, 1)  and  each  R  >  0  the  set  U  Fr  R  is  precoopact  in 

n 

*  *  hJ(T). 

The  concluding  result  of  this  section,  in  which  the  H1  topology  still  plays  the 
major  role,  is 

Theorem  5.9.  Let  5  e  (0, 1)  and  R  >  1  be  fixed.  Then 

(a)  The  set 

®«,R  "  lia  Fn,«,R 
n 

is  a  compact,  connected  subset  of  B(R)  (defined  in  4.4)  with  (1^,0)  e  G^R. 

(b)  Bach  (X,w)  in  G^  R  satisfies  parts  (a),  (c),  and  (d)  of  Theorem  5.1  and 
Vw(x1,x2)  t  Qj  for  all  (xvx2)  e  T-  C  T+. 

(c)  If  (X,w)  t  Gj  R\{(Xd,0)},  then  X  <  Xd,  w  >  0  in  T,  and  3w/3x1  <  0  on 
(0,«)  x  I. 

Proof,  (a)  Froa  Leama  4.1  and  Theorea  5.8  one  can  conclude  that  G^R  is  a 
compact,  connected  set  in  B(R).  It  was  shown  in  Lemma  5.5(c)  that  un  ♦  as  n  ♦  ■ 

and  thus  (Xd,0)  e  G^  R. 

Parts  (b)  and  (c)  follow  froa  Theorems  5.1,  5.4,  5.6,  and  5.7. 


q.e.d 


6.  The  Set  6 


j  R  as  a  Subset  of  *  *  (hJ(T)  n  C0'  \t)) 


To  show  the  existence  of  the  unbounded,  connected  set  of  solutions  S  described  in 


pert  (e)  of  Theorem  2. 1  two  issues  must  be  confronted.  The  first  is  the  connectedness  of 

I,™ 


Gg  R  when  the  stronger  topology  of  hJ(T)  fl  C  '  (f)  is  used  in  place  of  the  Hq 

topology  for  w.  The  second  is  roughly  to  show  that  in  the  set  Gg  R  there  are  solutions 

for  which  w  has  size  R  in  Hq  or  size  6  1  in  C®'  \  If  this  can  be  done  for  any 

R  >  0  and  any  5  e  (0,1),  then  the  unboundedness  of  S  will  follow.  It  will  be  shown 

that  in  G.  _  there  is  a  pair  (X,w)  for  which  |w|  .  -  R  or  for  which 

0  ,K  u  1 

Vw(x1#x2)  £  3Qg  for  some  (x^Xj).  This  almost  resolves  the  second  issue.  The  final 
step  is  to  show  that  for  sufficiently  small  6 ,  Vw  cannot  take  values  on  that  part  of 
3Qfi  where  p2  =*  “1  ♦  6. 


Theorem  6. 1.  For  each  S  e  (0, 1)  and  each  R  >  0,  the  set  Gg  R  is  connected  in 
R  x  (Hq(T)  D  C0'  \t)>. 

Proof.  Since  Gg  R  is  connected  in  the  original  topology  of  R  x  hJ(T),  it  will 
suffice  to  show  that  on  Gg  R  the  stronger  topology  coincides  with  the  original  one. 

That  is,  if  solutions  (Xk,wfc)  of  the  problem  with  discontinuous  density  pd  converge  to 

.  o,1- 

(X,w)  in  R  x  Hq,  then  w^  converges  to  w  in  C  (T).  The  HBlder  estimates  in 

Theorem  5. 1  are  satisfied  by  each  wR,  uniformly  in  k.  The  compactness  of  the 

0,  1  - 

collection  {wR}  in  C  (T)  is  then  assured  and  the  desired  convergence  follows 
immediately.  q.e.d. 


Recall  that  for  the  density  pn  the  set  of  solutions  n  emanating  from  (un,0)  is 
unbounded  in  R  x  hJ(T)  and  that  Fn# j  R  1«  the  maximal  connected  subset  of 
{(X,w)  e  En|Vw(x1fx2)  e  Qfi,  (xvx2)  e 

which  contains  (u  ,0)  and  is  contained  in  B(R).  He  can  assume  n  is  large  enough  so 
n 

that  y  <  2X  and  hence  for  each  n  there  is  a  pair  (X  ,w )  e  F  ,  _  for  which  either 
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*) 


B)  7wn(x<1n),x‘n))  €  8Q4  for  SOM  (x(1n>,x|n)>  e  T  . 

It  alternative  A)  hold*  for  an  infinite  eequence  of  valuee  for  n,  a  subsequence/ 

still  denoted  by  (^n'wn)  may  be  chosen  so  that  Xn  ♦  X  and  wn  converges  weakly  to 

w  in  hJ(T)  as  n  ♦  •.  Theorems  5.4  and  5.7  show  the  convergence  to  be  strong  in  hJ 

and  hence  there  is  a  pair  (X(w)  in  G.  _  with  Iwl  .  -  R. 

5»"  Hq 

If  alternative  B)  occurs  infinitely  often  it  can  be  assumed  that  a  subsequence  is 
chosen  so  that  the  following  result  is  applicable. 


Theorem  6.2.  Suppose  (X  .w  )  e  F  .  n  *  1,2,...  and  that  for  each  n  there  exists  a 
■  1  n  n  n/O/K 

point  (x^n*,x^n*)  in  T  for  which 

Vwn(x(1n),x£n)>  e  9QS  .  (6.1) 

Then  there  is  a  pair  (X,w)  e  G^  R  and  a  point  (x.j,x2)  such  that 

TwCx^Xj)  e  9QS  ,  (6.2) 

where,  if  x^  “  0,  the  limit  from  above  or  below  *>  0  is  intended. 

Proof.  The  proof  is  divided  into  three  cases  according  to  which  part  of  the  boundary 
3Q4  is  in  question.  Once  again  a  subsequence  can  be  chosen  and  renumbered  so  that  wR 
converges  to  w  in  Bq(T),  Vw  c  Oj  on  T,  and  one  of  the  following  three  cases  is 


relevant. 


easel.  ||i-  (x<"’,x‘n,)| 


5  t  n  -  1,2,  i 


First  not*  that  Corollary  4# 7  qua ran tees  that  Vwft  decays  to  sero  at  infinity# 
uniformly  in  n,  so  that  (x<1n),xjn>)  remains  bounded.  Without  loss  of  generality 
assume  that  these  points  converge  to  (x^,x2>.  From  the  HSlder  estimate  (4.20)  it  is 
clear  that  Sw^/Sx  1  converges  uniformly  to  3w/3x  ^  and  thus 


||2-  (S,,x2)|  »  «_1 


implying  (6.2),  since  Vw  t 


•  e.T  •  *  .  *  -  '  **;  •*,  »*,  . 


Case  2.  (xyw,x$“')  -  «  'in  -  1,2, 


As  with  Case  1  it  can  be  assusMd  that  (x^n*»x£n*)  ♦  (x^x^).  Zf  x2  f  0, 


the  H81der  estimate  (4.21)  applies  to  give 


■I—  (x  ,x  )  “  6  1 
3x2  1  2 


and  hence  (6.2). 


Suppose  Xj  "  0.  From  Lenta  4.3  it  can  be  aasuated  that  both  wR  and  3wn/3x.j 
converge  uniformly  in  a  neighborhood  N  of  (x^,0).  The  functions  3wn/3x2  cannot,  in 
general,  converge  uniformly  since  3w/3x2  is  discontinuous  along  (x2  “  0} .  However, 
from  (5.3),  the  limits  of  the  x2  derivatives  of  w  from  above  and  below  {x2  “  0) 


exist  and  we  define 


<2  “  V<;1'X2> 

2  x2+0+  i 


w  (x,,x.) 


*2  x_>0-  *2 


The  link  among  these  various  limits  is  provided  by  the  combination 

Pn(x2)(f2(Vwn)  -  Xnwn) 

which,  according  to  Theorem  4.4b),  has  a  HSlder  exponent  a  >  0,  uniformly  in  n.  This 


information  combines  to  give 


11m  p  (x‘n,){f,(wir  ,«_1)  -  Xw) 
n  n  2  2  x, 


5  -  Xw> 

*  X1  x2 


-  P_{f,(w  ,w  )  -  Xw) 
X1  x2 


where  the  terms  other  than  w_  are  evaluated  at  (x,,0). 


Now  one  considers  cases.  If 


f,(w  ,«  )  -  Xw  -  0  , 

x, 


'  V  >  >  •  »  •  *  "  ■  *  •  »  */  v  J  <T  W.  ^  *  V  N; 


t 


.■»  j  wsrw^rwBmmasKsarrcvwmv3&WB mg  u.n  n. 


■W- 


than  the  bracketed  terms  in  (6.4)  and  (6.5)  must  also  vanish,  yielding 

f2<w*  <*C  )  “  )  “  f2<Wv  '*  ^  * 

*  *1  «2  ^x1x2  2x1 


Since  f 22  >  °»  it  follows  that 

and  thus  (6.2)  holds* 

Next/  suppose 


»“  •  w*  -  4  1  » 
x2  x2 


f,(w  ,6-1)  -  Xw  >  0  . 
*  X1 


(6.6) 


Since  p.  <  p  (x,)  for  all  x,,  it  followa  from  (6.3)  and  (6.4)  that 
▼  n  *  < 


and  so 


P+*f2("*  '4-,)  "  *  P+{f2<W*  **C  1  "  Xw* 

T  *  T  4  *2 


f,(w  *8-1)  <  f2(w.  »w*  >  • 

2  X  ^  2X^X2 


Again  using  f22  >  0  ons  concludes  that 


<  w! 


and  ainca  any  limiting  values  of  the  gradient  of  w  must  lie  in  Qg,  (6.2)  follows. 

If  the  inequality  in  (6.6)  is  reversed  an  argument  similar  to  that  juat  given,  but 
uaing  the  equality  of  (6.3)  and  (6.5),  leads  to  (6.2). 

.,2 . 

^  -  1 
The  argument  here  is  similar  to  that  given  in  Case  2,  but  with  6  replaced  by 

•1  +  6.  q.e.d. 


Our  next  task  is  to  show  that  w_  cannot  assume  the  value  -1+6  if  6  is 

*2 

sufficiently  small.  Sines  -1  +  6  <  w  <  6  \  the  value  -1  +  6  would  be  a  minimum 

*2 

for  w_  . 
x2 
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is  assumed  on  ths  lins  {xj  «  0}  as  a  limit  from  above 


T  emma  6.3.  The  minimum  of  w„ 


or  below. 


Proof .  Since  w  satisfies  the  equation 


•j|-  f^Vw)  -  0  in  T+  U  T" 


it  follows  from  differentiation  that 


3  .  t  n  \  3  3w  ,  _+  _■ 

- —  f  .  .  ( Vw)  - —  - —  -  0  in  T  U  T 
3xa  ij  3Xj  Sx^ 


for  k  “  1  or  2.  By  the  maximum  principle,  w  ,  which  is  continuous  in  T  and  in 

x2 

T~,  must  take  its  minimum  value  on  the  boundary  of  T+  U  T~.  Since  w  is  positive  in 

T  and  vanishes  for  x2  ”  -h  and  for  x2  “  1  -  h,  the  strong  maximum  principle  implies 

that  wx^  >  0  for  x2  “  -h  and  wx^  <  0  for  x2  “  1  -  h.  Consequently  the  minimum  is 

negative  and  must  occur  where  x2  •  0  or  x2  »  1  -  h. 

The  function  w_  is  negative  in  (0,«)  x  I+  (cf.  Theorem  5.4)  and  is  sero  on  the 
X1 

half  line  (0,»)  x  {1  -  h}.  Since  w  is  even  in  x.,  w  is  sero  on  {0)  x  I+  and  the 

*1 

etrong  maximum  principle  implies 

w„  >0  on  (0,-)  x  {1  -  h}  (6.9) 

*1*2 


w  <0  on  {0}  x  I+  . 
*1*1 


(6. 10) 


The  inequality  (6.9)  implies  that  the  minimum  of  wx^  on  the  line  x2  “  1  -  h  must  occur 
at  x.  ”  0.  Suppose  the  equation  (6.7)  is  written  out  in  nondivergence  form.  Since 


w  «  0  on  {0}  x  I+,  it  follows  from  (6.10)  and  the  ellipticity  of  (6.7)  that 
*1*2 


w„  _  >0  on  {0}  x  I+  . 

x2*2 


(6.  11) 


Hence  the  minimum  of  w  on  {0}  x  i+  must  occur  at  (0,0),  showing  that  the  minimum 

x2 

must  occur  on  {x2  “  0}. 


y . \  ■■  v*y  'A  ■  \ A'  .’•yA  /■  j  ••  v-v  n  ‘."-V"  >■- 


The  ntxt  result  establishes  a  lower  bound  for  X  which  aervaa  in  Theorem  2. 1  part 

(b)  and  which  will  also  be  used  in  obtaining  a  lower  bound  for  w„  .  The  proof  of  thie 

x2 

lemma  will  be  given  in  connection  with  related  computations  at  the  end  of  Section  8. 


Leasts  6.4.  For  any  6  c  (0,1)  and  R  >  0  let  (X,w)  belong  to  G^  R.  Then 


X  > 


1  (^P_  + 

2"  P_  ~  P  + 


(6.  12) 


Recall  that  the  awitch  from  spatial  to  aemi-Lagrangian  independent  coordinatea  waa 
Bade  to  circumvent  the  difficultiea  Inherent  in  a  direct  approach  to  a  free  boundary.  For 
each  aolution  (X,w)  e  Gj  R  the  free  boundary  ia  merely  the  graph  r  •  {x,w(x,0)>  of  the 
function  w(x,0).  Recall  from  Theorem  5.9  that  w(«,0)  c  C1,a(R).  in  Section  7  it  will 
be  shown  that  it  is  real-analytic.  Having  obtained  the  boundary  separating  the  two  fluid 
components  we  now  revert  to  the  paeudo-atream-function  f(x,y)  where  (x,y)  are  spatial 
coordinatea  in  the  strip  S  (cf.  Section  1).  Recall  from  (3.10)  that 


in  S+,  the  upper  fluid  region,  and  in  S~,  the  lower  one,  reapectively .  Since 
c2  >  g/Xd»  (6«  13)  will  give  a  lower  bound  for  wx^  if  an  upper  bound  for  ia 

establiahed.  The  next  lemaut  is  a  first  step  towards  bounding  . 


Lenma  6. 5.  Let  (X,w)  belong  to  G^  R  and  let  be  the  corresponding  pseudo-stream- 
function.  Then  for  any  real  m 


m+1  , 

/  /  |V*|  dxdy  <  c2!^  (6.14) 

m-1  I 


where  depends  only  on  and  h. 


Remark.  Froa  T  >—■  6.4  it  la  clear  that  c^K ^  -  gK^/X  la  bounded  in  terms  of 
and  h. 

Proof  of  Leaaa  6.5.  Let  the  integral  in  (6.  14)  be  denoted  by  E.  From  the  relatione 
(2.20),  (2.30),  (3.  10),  and  (3.11)  there  reeulta 

■♦1  ?( 1-h)  1  +  y2 
«  -  /  /  — s-5  <*<* 

■-1  'T(-h)  y* 


.  ,_v.  1  +  w 

m+ 1  1-h  x . 

*/  /  T+"*  c/<5d(X2)  **1  dxT  **2 

i-1  -h  x2  2 


m+1  1-h  1  +  wx, 

“=2/  ,  /  OdwV 

■-1  -h  x. 


dx jdx2 


(6.  15) 


Let  u(x.j)  be  a  saooth  cutoff  function  which  equala  1  on  [■  -  1,a  ♦  1]  and  haa 

support  on  [a  -  2, a  +  2).  Consider  the  fundaaental  equation  (4.1)  for  a  pair 

2 

(>n,wn)  e  Fn<jfR  converging  to  (X,w)  c  G5R.  Multiplying  (4.1)  by  u  »n,  integrating 
by  parts,  and  letting  n  ♦  •,  one  obtains  (cf.  4.11) 


J  ff 


2  ♦  w 


2  1  '  d  v  2 

*  t  (i  +  w  r 
“2 


2  )|Vw|2 


X1  “22 
-2  //  0duu'  f  +  w — ♦  X(p_  -  P+)  /  U  (x, )w  (x1 , 0 )dx 


2  *  1 

<  c  //  p  U  - -  + 

T  ( 1  +  w  ) 

*2 


P_  *+2  , 

,  —  /  ( W *  >  dx  +  4X(p_  -  p  ) 

e  a-2 


(6. 16) 


where  |w|  <  1  (cf.  4.35)  has  been  used.  If  e  is  set  equal  to  1/2  in  (6.16)  the 
resulting  inequality  implies 
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2  ‘  I  pd  . - T  X  T  2  >  j  pd  1  +  w  ' 

B-1I  ( 1  +  w  )  2  m-1I  x, 

X2 

where  K<  depends  only  on  and  h.  The  last  Integrand  in  (6.  15)  contains  a  sun  of 

(1  +  w„  )~ 1  and  w?  (1  +  w„  )*'.  The  second  of  these  also  occurs  in  the  second  integral 
x2  x  1  *2 

of  (6.17)  and  so  is  estimable. 

Now  consider  p.(1  +  w„  )-1.  Where  w„  >  -  4,  this  integrand  is  at  most  2p  . 

a  *2  x2  *  a 

while  on  the  complementary  set 


2  +  w 


- - -  <  i 

1  +  w_  3 


(1+w  ) 


2  2 

- w 

2  x2 


Use  of  these  estimates  with  the  first  integral  in  (6.17)  gives  a  bound  for  the  integral  of 


p.(  1  +  w  )” 1  and  the  combined  estimates  yield  an  estimate  for  (6.15).  q.e.d. 


To  make  use  of  the  bound  (6. 14)  it  will  be  necessary  to  show  that  the  distance  from 
the  streamline  T  to  the  boundary  88,  denoted  dist(r,3S),  is  bounded  below. 


Lemma  6.6.  Suppose  (X,w)  e  Gj  R.  Then 

dist(r,8S)  >  K  >  0  (6.  18) 

where  K  depends  only  on  p^  and  h. 

Proof.  In  the  case  of  elevation  waves  (e  >  0)  being  treated,  w  >  0  so  the 
distance  from  T  to  the  lower  boundary  is  at  least  h.  Since  w  is  even  in  x  i  and 
decreasing  on  0  <  x^  <  •»  the  distance  to  the  upper  boundary  is  d  •  1  -  h  -  w(0,0). 

Consider  polar  coordinates  centered  at  (x,y)  “  (0, 1  -  h)  so  that  points  within  the 
strip  S  have  angular  coordinates  in  (0,w).  From  Lemma  6.5 

lx  -  1-h  1  2 

/  /  | V*|  rdrdB  <  /  /  | 7*1  dxdy  <  C2K, 


(6.  19) 


where  u+  “  sup(u,0)  and  u~  »  inf(u,0).  then  f(r)  is  an  increasing  function  of  r 


for  0  <  r  <  R. 


Theorem  6.8.  Suppose  (X,w)  e  G^  R.  W»en 


(6.24) 


(6.25) 


where  the  constants  depend  only  on  and  h. 


Proof.  In  Lemma  6.3  it  was  shown  that  the  minimum  of  1  +  w_  is  assumed  on  the 
-  x2 


line  {x2  ”  0),  from  above  or  below.  From  the  discussion  surrounding  (6.13)  an  upper 


bound  for  3t|i/3y  on  the  upper  and  lower  sides  of  r  will  yield  (6.24).  The  estimate 


(6.25)  follows  similarly  from  e  bound  on  3<i/3x  which  is  harmonic  in  S  and  must  assume 


its  maximum  and  minimum  on  3S+  U  3S  .  Since  la  constant  on  the  lines  (y  m  ”h)  and 


{y  *  1  -  h},  it  suffices  to  bound  |<i  |  on  T. 


Let  i|>+  denote  the  restriction  of  4*  to  S+  and  ,  the  restriction  to  S“ 


Since  i|»+  >  0  in  S+  and  $  <  0  in  S”,  this  notation  is  consistent  with  that  in  the 


previous  lemma.  According  to  Theorem  5.1  the  function  w(x,0)  describing  T  is  of  class 


C1'0  and  Vw  is  of  class  C°  in  T+.  Since  Vw  lies  in  ,  the  relations  (3.10)  and 


(3.11)  show  that  V*4  is  of  class  Ca  in  S*,  respectively.  Let  (x,y)  be  a  point  on 


T  and  let  c*  "  {v^*(x,y)|.  If  sr  denotes  the  ball  of  radius  r  centered  at  (x,y>. 


then  from  the  previous  lemma 


(r)  -  -j  /  +  |V*+|2  •  -y  /  .  |V* 


(6.26) 


r  B  ns 
r 


r  B  ns 
r 


k>.v;>.v.v.vvV  ■  - 


.  •.  ■-  *-  •-  .  •-  •  c  v  > .  •h.  >  .  v'V  . -  o'  - 

,**  •  •  *  .  ^  "  e  I \  •  •*  •  M  ■’Jm  “  J,  m  W  w  — %  L  V  W  w 


is  increasing  in  r  for  0  <  r  <  R  where  R  »  min(h,  dist(r,3S)).  From  the  regularity 
just  cited  it  follows  that  f(r)  has  a  limit ,  denoted  *(0),  as  r  ♦  0  and 

~  (a*)2  •  \  (a~)2  -  f(0)  <  f (R)  <  (R^c2*,)2  (6.27) 

by  Lensa  6.S. 

Recall  that  the  continuity  of  pd<f2(7w)  -  Xw)  (cf.  4.31)  is  a  translation  of  the 
pressure  condition  (2.28)  which  becomes 

j  lie*)2  "  (c  )2]  -  (p  +  "  P _ ) t— j  -  gy>  (6.28) 

in  the  current  notation.  The  two  relations  (6.27)  and  (6.28)  easily  imply  that  a*  and 
o  are  bounded  in  terms  of  p^  and  h.  q.e.d. 


It  follows  from  the  last  result  that  there  is  a  SQ  >  0  so  that  for  0  <  6  <  SQ  a 


R  >  0  there  is  a  pair  (X,w)  e  G  for  which 

0  0  R 


1*1  ,  -  R 

“0<T> 


«ax( J w 


*l»x  I  )  -  «_1  • 

1  m  2  • 

L  (T)  L  (T) 


Hence  S  C  R  *  (hJ(T)  n  C°'1(T))  defined  by 

S-UG, 

n  —  ,n 
n 


is  an  unbounded,  connected  set  of  solutions  of  problem  ( 2. 36)- ( 2* 39) 


7.  Proof  of  Theorem  2. 1 


7.  1.  Existence  end  Besic  Properties 

The  major  part  of  the  work  involved  in  proving  Theorem  2. 1  has  been  done.  In  this 
section  we  summarize  those  results  already  obtained  and  add  proofs  for  the  remaining  parts 
of  the  theorem,  with  one  exception.  As  already  noted,  the  lower  bound  X_  for  X,  given 
by  (6.12),  is  derived  in  Section  S. 

Part  (a).  The  basic  existence  result  in  Theorem  2.1  was  the  culmination  of  the 

analysis  in  Sections  3  through  6  and  waa  asserted  at  the  end  of  the  previous  section. 

That  elements  of  S  are  solutions  of  (2. 36)-(2. 39)  follows  from  Theorem  5.1(d).  The 

inclusion  of  (X.,0)  in  S  follows  from  Lemma  5.5(c). 
a 

Part  (b)  The  inequality  X  <  X^  was  given  in  Theorem  5. 1.  From  Theorem  5.6  it 

follows  that  (X.,w)  e  5  implies  w  =  0.  All  functions  w  under  consideration  have  been 

d 

even  in  x 1  and  hence  their  limits  are  even  in  x^.  The  remaining  parts  of  b)  were  shown 
in  Theorem  5.4. 

Part  (c)  The  analyticity  of  w  in  the  closed  strips  T*  will  be  shown  in  Theorem 
7.6  below.  As  for  the  other  assertions  of  regularity,  they  follow  from  Theorem  5.  1. 

Part  (d)  These  bounds,  which  tranalate  into  bounds  on  velocities,  were  obtained  in 

Theorem  6.8. 

Part  (a)  In  Im«sa  7.4  of  [26]  the  exponential  decay  of  solutions  was  shown  both  for 
periodic  waves  over  a  half  period  and  for  solitary  waves.  The  proof  in  [26]  can  be 
simplified  if  only  solitary  waves  are  considered. 

Part  (f)  This  nonexistence  result  is  shown  in  Theorem  7.7  below. 

In  the  event  e  <  0,  one  can  use  (5.43)  to  arrive  at  the  reversal  of  inequality 
(5.47)  and  so  prove  Theorem  5.6.  Similar  small  changes  are  required  for  the  proof  of 
Theorem  5.7. 
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f. 


7.2.  Analyticity  of  the  Interface 

As  was  observed  by  Klnderlehrer,  Nirenberg  and  Spruck  [ 15]  the  analyticity  of  w  in 
T*  will  follow  from  results  of  Morrey  ([20];  [21],  Chap.  6)  once  it  is  shown  that  w  is 
of  class  C2  in  T*.  In  [  15]  solutions  of  free  boundary  problem  are  examined  under  a 

2 

typical  assumption  that  the  solution  is  known  to  be  C  .  Then  a  partial  hodograph  or 
Legendre  transformation  together  with  suitable  combinations  and  reflections  of  solutions 
in  abutting  domains  produce  an  elliptic  system  adapted  to  Morrey' a  results.  Here  we  have 

worked  in  a  semi-Lagrangian  (partial  hodograph)  setting  for  most  of  the  analysis  and  hence 

2  +  2 
need  only  do  a  reflection  once  the  C  character  of  w  in  T  is  established.  The  C 

estimates  could  be  carried  out  using  the  techniques  from  [3 ] ,  but  a  relatively  short 

direct  proof  can  be  given  and  we  do  that  in  Corollary  7.5.  A  lemma  is  required  (and  could 

be  given  in  a  local  version). 


Lemma  7.  1.  Let  (A,w)  be  an  element  of  Gj  R  (cf.  4.7).  Then  there  is  a  p  «  p(5  )  >  2 
and  a  constant  K  depending  on  6  and  R  such  that 


/  v«  <  K{1  +  (/  | Vv  |2,2-<P/2>j 


T  1 


for  i  -  1,2,  where  v  -  w  . 

*1 

Proof.  It  will  suffice  to  work  with  an  element  (A,w)  e  F  .  „  (see  Section  4)  and 
—  n,o  ,k 

obtain  an  estimate  which  is  independent  of  n.  The  subscript  n  will  be  suppressed. 

The  function  v  satisfies 

■  t> r  (pfii<7w’  f£->  *  "Xp'v  ln  T  (7,2) 

and 

v  -  0  on  3T  .  (7.3) 

From  inequalities  (4.9)  and  (4.29)  it  follows  that  there  is  a  p  -  p(fi)  >  2  such  that 

| Vv|  _  <  k|w|  ,  <  K'R  (7.4) 

LB(T)  L  (T> 

for  all  s  ln  the  interval  [2,p(6)J,  where  K  and  K'  depend  on  5 .  A  version  of  the 


>vv 


Gagliardo-Nirenberg  Inequality  will  be  needed  for  a  function  u  e  H'(T).  By  extending 

u  as  an  even  function  across  the  lines  {x2  -  -h}  and  {x2  -  1  -  h>  and  using  a  cutoff 

function  in  the  x2  direction  one  easily  derives 

|»|  4  <  K1|u|1:(p/4,Mp/4  (7.5) 

L4(T)  11  H’tT)  lP(T) 


from  the  main  theorem  in  lecture  II  of  [22].  Here  p  can  take  any  value  in  the  interval 
[1,4].  He  can,  without  loss  of  generality,  assume  p(6)  <  4  and  let  p  -  p(6). 

Equation  (7.2)  can  be  written  as 

sr-  !r-  +  If-  -  *pv>  -  -  s=-  <p* „  -  xp  |f- 


3x2  21  8x,  K  22  8x2 


8x.  1j  8x. 


fijk  denotes  a  third  derivative,  and  the  summation  convention  is  used.  Without  using  the 
equation  one  has 

8x^  <pC2i  aTj-  +  Pf?2  8^  "  Xpv> 


-  pl‘21k  fe  fc  +  <21  ^ *'22k  ^  **22  1^  - 


With  pf_.v  -  Xpv  denoted  by  u,  it  follows  from  (7.6)  and  (7.8)  that 

*3  * j 

f  1 2  -  r  r  r  / 1  r?..  1 2  .  I  o _ I  2  .  4  .  n 


/  |Vu|2  <  Kjt/  <|vYr  ♦  \Vv\  *  v*  ♦  v*  )] 

mm  12 


where  K2  depends  on  i  and  R.  Since  the  L2  norm  of  | 7v|  and,  by  embedding,  any 
if  norm  of  v  are  bounded  in  terms  of  R  (cf.  (4.9))  such  terms  can  be  accommodated  by 
including  additive  or  multiplicative  constants.  Since  f 22  is  bounded  below  by  a 
positive  constant  depending  on  5,  the  use  of  (7.5)  yields 
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T  2  T 


*4  /  «4  +  *',</  y4  +  /  "J  > 


T  T  1 


<  |u|p  +  1  +  /  v*  ] 

l  (T)  LP(T)  T  X1 


(7.  10) 


Using  (7.4)  to  estimate  u  in  lP  and  (7.9)  to  estimate  Vu  in  L2,  we  have 

/  \  <  V1  *  f  \  +  (/  <1  vyI2  +  v*  +  v*  )>2'{p/2,i  . 

T  *2  T  X1  T  X1  X2 


(7.  11) 


Since  2  -  (p/2)  <  1,  equation  (7.11)  provides  a  bound 


/  \  <  vi +  f  *  {f  I*yI2>2"<p/2>] 


T  1  T 


(7. 12) 


where  K?  depends  only  on  6  and  R. 

The  embedding  inequality  (7.5)  applied  to  v„  e  hJ(T)  gives 

X1  u 


f  v*  <  *B<f  |vr|2>2"<p/2)  /  |v  |p  <  K  {/  |Vy|2}2"(p/2) 

T  1  T  T  X 1  *  T 


(7.13) 


again  using  p  -  p(5)  and  (7.4).  The  use  of  (7.13)  alone,  and  in  (7.12)  produces  the 
result  (7.1).  q.e.d. 


Theorem  7.2.  Suppose  (X,w)  <  G.  .  Then 

0|R 


/  |7w  I  <  K 
T  X1X1 


(7. 14) 


where  K  depends  on  2  and  R. 


Proof.  Continuing  the  analysis  in  the  previous  proof  for  a  solution  in  F  ,  _  one 

n,o ,  R 


differentiates  (7.2)  to  obtain 


-^7  <pfij(7w)  §->  (pfijk  "  Xp’Y  ln  T 


ij  ijk  3Xj  9x)t 


(7. 15) 


a  .v/Av.sv-  w>v- 


yy  *  v  /  «* 


2  2 

and,  since  y  ■  3  w/dx^ 


Y  “  0  on  8T 


The  ellipticity  of  the  equation  (7.15)  yields 

v  I  |7y|2  *  L  pfij  It  It 

T  T  i  j 


“  "  l p  f ijk  %  +  2X  l  PTY» 


<  y  /  i^yI*  ♦  ;  <vj  ♦  vj  ♦  y2) 

T  T  1  2 


(7.  16) 


(7.17) 


The  bound  (7.14)  now  follows  frost  ths  previous  lesata  and  T asms  4.2.  We  note  again  that 


the  estimates  are  all  uniform  in  n  and  so  hold  in  G. 


q.e.d. 


Corollary  7.3.  If  (X,w)  e  G^  R,  than 


(7. 18) 


1  LP(T) 


for  all  p  e  [2,»),  where  K  depends  only  on  p,  6,  and  R. 

Proof.  Proa  the  previous  theorem  the  integral  of  |  Vy  | 2  over  T  is  bounded  and 

from  Lemma  7. 1,  so  are  the  integrals  of  vj  and  v*  .  Proa  (7.9) 

*1  *2 


/  l7<pf21-5t+  pf22  it  "  *PV)|  <  K1 

T  12 


where  K1  depends  upon  6  and  R.  Hence  P£2ivx  “  *PV  ia  in  iP (T)  for  all 
P  e  (2,«),  according  to  embedding  theory.  The  same  embedding  result,  used  with  (4.9)  and 
Theorem  7.2,  shows  that  v  and  vx^  are  in  lP(T)  for  all  p  e  (2,») .  Since  f22  Has 
a  positive  lower  bound  depending  only  on  6,  it  follows  that  v  is  in  the  saae  I P 


spaces. 


q.e.d. 


It  is  now  possible  to  return  to  the  elliptic  equation  for  y  to  obtain 


y*  .*•  *  .**  “•*  *.'■  "•  eV*  *•'  "«  N*  ‘/l  V  •  v  •' * 


i^eeest 


Theorem  7.4.  Let  (X,w)  be  in  G,  _.  Then  there  is  a  q  “  q(6)  >  2  so  that  y  m  w 


satisfies 


(7.  19) 


■  w.  _  can  be  written 


/  |Vy|q  <  K 

T 


where  K  depends  on  fi  and  R. 

Proof.  For  (X,w)  e  F  .  _  the  equation  (7.15)  for  y 
1  n, o ,k 

-  ^  "'u  '  * h 


F1  pC1jk  3Xj  3^  ' 

P2  “  pf  2jk  li“  fe-  “  XPY 


h-Xpjl- 

a*2 


By  the  previous  corollary  F,  and  F2  are  in  tP(T)  for  p  e  [2,»)  while  from  Theorem 
7.2,  he  L2(T).  As  for  equation  (4.24)  Theorem  1  of  Meyers  [18]  yields  (7.19)  for  some 
q  >  2.  All  estimates  are  independent  of  n  so  the  inequality  (7.  19)  follows  for  all 


(X,w)  e  G, 


q.e.d. 


2  a  ^ 

Corollary  7.5.  Let  (X,w)  be  in  G^  R.  Then  w  e  C  '  (T  )  for  some  a  >  0  and  hence 

2fo 

the  fluid  interface  y  “  w(x,0)  is  of  class  C 
2  ci 

Proof ■  The  C  '  character  of  the  interface  follows  from  (7.19)  and  embedding 
theory.  The  Schauder  theory  ([13],  Lemma  6.18)  provides  the  same  regularity  in 


Theorem  7.6.  Suppose  (A.w)  is  a  solution  of  (2. 36)-(2. 39)  belonging  to  S.  Then  w  is 
real  analytic  in  T~  and  T^. 

2fa  + 

Proof.  Since  w  is  of  class  C  in  T* ,  the  techniques  of  [15]  can  now  be 
used.  In  fact,  if  we  revert  to  the  formulation  in  which  f(x,y)  is  harmonic  in  S+ 
and  S”  and  satisfies  (2.27).  (2.28).  then  nieorem  3.2  of  [15]  yields  analytlcity. 

7.3.  Nonexistence  of  Shall  Solutions  with  e  •  w  <  0 
We  now  prove  part  (f )  of  Theorem  2. 1. 

Theorem  7.7.  Suppose 


e 


>  0 


There  is  a  positive  n  such  that  if  (X.w)  c  1  x  (hJ(T)  n  C0' *(T))  is  a  solution  of 
(2.39)  with  w  <  0,  |2-  >  0  for  x,  >  0,  and 


|l  -  A. |  +  | w |  .  + 

d  hJ(t) 


|w|  Q  *  <  n 

C°''(T) 


(7.20) 


then  w  =  0. 

Proof.  Repeating  estimates  already  used,  in  conjunction  with  difference  quotients  in 
the  x  i  direction,  one  can  show  w  is  smooth  in  each  of  T~  and  T^  and  we  shall 
assume  the  smoothness  necessary  for  the  calculatione  to  follow.  First  it  will  be  shown 
that  if  the  hypotheses  of  the  theorem  are  met  and  X  >  X^,  then  w  2  0.  The  steps 
involved  in  the  proof  of  Theorem  5.6  up  to  (5.40)  can  be  duplicated  for  the  case  of  a 
discontinuous  density,  that  is,  for  p^.  The  equation  (5.  12)  which  formally  becomes 

<PdY')'  -  Xp^y  -  opdY 

is  replaced  by  the  weak  form 

-  /  p  .y'f '  +  X  /  Pd(Y?)'  “  a  /  PdY*  (7.21) 

I  I  I 
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vher*  y  f  are  in  Hq(I).  Equation  (7.21)  arises  from  the  variational  problem 

tX(p_  -  p+)u2(0)  -  /  Pd(u*)2} 

a  »  max  - - 

u«hJ(I)  /  Pd« 

which  correspond,  to  (5*11).  The  maximizer  y  is  normalise  by  /  Pdf2  ”  1  and,  as  in 

the  proof  of  T easts  5.5,  we  have  a  >  0  since  X  >  A  ..  In  following  the  proof  of  Theorem 

<1 

5.6  one  need  only  let  the  derivative  of  a  density  in  an  integral  be  interpreted  as 

/  p)Jg(x1,x2)dx2  -  -(p_  -  p+)g(x1,0) 

in  the  case  of  the  discontinuous  pd  from  (2.26).  As  regards  (5.35),  let  t  correspond 
to  the  density  pd*  The  extra  condition  of  orthogonality  yields  r  >  Ad  and  it  can  be 

assuawd  that  X  is  sufficiently  close  to  X.  to  make  1  -  X/t  ■  8  >  0. 

As  regards  the  steps  in  the  proof  of  Theorem  5.6  past  (5.40)  the  only  change  is  in 

the  discussion  following  (5.40),  for  now  a  is  positive  rather  than  negative.  Inequality 
(5.41)  is  replaced  by 

<f')2  +  /  pJvrI2  <  ^f-r2  +f  o*  ♦  k-f3 

and  in  place  of  (5.42)  one  has 

G'  +  or  ♦  er2  >  -2eF2  -  f2  +  u*  k*f3)  . 

It  follows  from  this  inequality  that 

G*  -  KU*  ♦  («  -  K.|f|  )r  >  0 

l" 

for  constants  K  and  (we  may  assume  |f|  m  <  1).  If  this  last  inequality  is 

1> 

integrated  from  x  to  where  x  >  0,  and  the  result  rewritten  in  the  vein  of  (5.45) 

the  result  is 

*x 

-  I  V  TTT  (x,x.)  -  T}dx2  +  <«  -  i  >  0  • 

X  X.  2  L  X 
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If  w  la  aaall  in  C0' 1 (T) ,  tha  expression  1  -  Rs/y  will  be  poaitiva  aa  will 
a  -  K,|r|  .  since  -w  and  F  ara  both  negative,  a  contradiction  ariaaa  unleaa 

1  j,  X1 

w  3  0. 

How  wa  nay  aaauae  X  <  Xd»  If  f  ia  aat  equal  to  w  in  tha  aquation  (2.39),  tha 
raault  ia 

|Vw|2(2+wx) 

/  Pd  - - y~  “  /  pdfi(7w)wx  “  X(p-  “  p+>  / 

T  a  2(1  +  w  )*  T  xi  — 

*2 

or,  rewritten, 

.  *x  <3  +  2-x  > 

/  pJVw|2  -  X(p_  -  p+)  /  w  (x1,0)dx1  -  /  pJVw|2  — - - .  (7.22) 

T  -•  T  2(  1  +  w  ) 

X2 

In  analogy  with  (5. 16)-(5. 18),  we  decompose  w  aa 

wtx^xj)  »  F(x1)Yd(x2)  ♦  R(xvx2)  (7.23) 

where  Y .  ia  the  function  deacribed  in  (5.13)  and 

Q 

F(x,)  -  /  p^U^XjJYjdXj,  i)tl. 

Mote  that  Yd  aatiafiea  the  weak  equation 

/  PflYdf '  -  Xd(p.  -  p+)Yd(0)»«>)  (7.24) 

for  any  f  e  Hq(i).  a  further  relation,  the  analogue  of  (5.19),  ia 

J  pdR(x j,x2)Yddx2  -0,  x1  e  R  •  (7.25) 


The  decomposition  (7.23)  and  the  relationa  (7.24)-(7.25),  uaed  in  the  left-hand  aide  of 
(7.22),  reault  in 
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/  / 


/  PdI Vw| 2  -  X(p_  -  p+)  /  W2(x1,0)dx1 


/  <F')  d*!  +  /  Pdl  VRI 2  -  X<P_  -  P  +  )  /  R2(x1,0)dx1 


E^v 

w 


2  2  m 
+  (P_  -  p+)Yd<0)(Ad  -  X)  /  F  dx.,  +  2(Ad  -  X)(p_  -  P+hd(0)  /  FR(x j,0)dx ^  .  (7.26) 


In  analogy  with  the  discussion  following  (5.35),  the  orthogonality  (7.25)  ensures  that 


V;.; 

fry' 


/  Pd!7Rl  -  x<p_  -  P+>  /  R2(x.,0)dx1  >  K  /  |vr|‘ 
T  —  1  T 


(7.27) 


for  sene  K  >  0  whenever  |x  -  Xd|  is  sufficiently  snail.  For  the  last  tern  in  (7.26) 
we  use  the  Schwarz  and  PoincarS  inequalities  to  obtain 


2(Xd  -  X)(p_  -  P+)7d(0)  /  F(Xl)R(x1f0)dx 


<  7  /  l7R|2  ♦  K*Ud  *  ^|2Td<°)  /  R2dx1  .  (7.28) 

T  -m 


If  (7.27)  and  (7.28)  are  used  in  (7.26)  and  if  Xfl  -  X  is  made  still  smaller  if 
necessary,  there  results 


i 


m 


Ul. 


I  PdlVwl  “  X<P_  -  P+)  /  w  (x1,0)dx1  >  /  (F'(X  ))2dx  +  j  f  | VR|  ‘ 
T  —  —  T 


Now  consider  the  right-hand  side  of  (7.22).  First,  for  any  e  >  0 


w  (3  +  2w  ) 

f  Pd"2  — - §-  <  e  /  |VR| 2  +  e  /  (F*  )2dx, 

Tax12(1+w)  T  — 


(7.29) 


(7.30) 


provided  w  is  sufficiently  snail  in  C0,1(T).  Next,  with  the  use  of  (7.23)  one  obtains 
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>  V\  V  s 


wx  (3  +  2wx  ) 

S  0A\  — - r  <  /  ed  *  I  <^(*1)Y,<x2»)3  +  K*  /  r4dx,  +  |  /  |7R|2  (7.31) 

T  2  2(  1  +  w  )  T  -*•  T 

x 
*2 

0,  1  - 

for  w  sufficiently  small  In  C  (T),  where  K  is  the  constant  appearing  in  (7.29). 

~  3 

The  first  term  on  the  right-hand  aide  of  (7.31)  is  the  integral  over  R  of  eF  (x1 ) 

where  e  is  defined  by  (5.26)  and  is  positive,  being  a  positive  multiple  of  e  in 

(2.41).  If  the  inequalities  (7.29).  (7.30)  and  (7.31)  are  used  with  equation  (7.22)  and 

the  integrand  F4  is  estimated  above  by  — J  F J  aF3,  there  results 

L 

_  ,  •  ,  "  , 

(J  -  e)  /  |7»r  +  (1  -  e)  /  (P-J-'dx,  -  (e  -  K* |r|  j  j  <  o  .  (7.32) 

T  -*•  L  -» 

Since  P  -  /  p.wy  .  <  0.  it  is  clear  from  the  estimates  just,  completed  that  if  w  has  a 
^  a  a 

sufficiently  small  norm  in  Hg(T)  nc0,1(T)  and  if  |x  -  X^l  is  sufficiently  email. 

P  =  |7R|  S  0  and  hence  w  5  0.  q.e.d. 


8.  Waves  of  large  Amplitude 

Since  the  branch  of  solutions  S  from  Theorem  2. 1  is  unbounded  in 
K  x  (Hq(T)  n  c0'  1  (T ) )  and  the  range  of  X  is  bounded,  the  norms  of  w  in  Hq(T)  and  in 
C0' ^  (T)  cannot  both  remain  bounded  on  S. 

In  Section  6  it  was  established  that  there  is  a  6~  >  0  so  that  w„  >  -I  +  j.  for 

U  *2  U 

all  (X,w)  c  S.  If  the  gradient  of  w  were  uniformly  bounded  for  (A,w)  c  S,  there 

would  be  a  positive  4 1  <  4g  such  that  for  all  solution  in  S,  V w  t  fij  .  But  then  the 

coefficients  f^j  would  satisfy  the  bounds  (4.15)  and  (4.16)  with  4  »  4  uniformly  for 

(X,w)  e  S.  This,  in  turn,  would  imply  that  in  2  there  would  be  solutions  having 

I w I  ,  arbitrarily  large.  For,  if  not,  then  the  local  estistates  (4.20)  and  (4.33), 
L2(T) 

2 

summed  over  m,  would  lead  to  a  bound  for  Vv  in  L  (T)  and  hence  to  a  bound  for  w 

1  2 
in  Hg(T).  TO  go  one  step  further,  since  0  <  w  <  1,  a  large  norm  for  w  in  L  can 

arise  only  by  having  a  wave  which  is  very  "broad”.  The  alternatives,  then,  are  to  have 

waves  which  are  arbitrarily  broad  or  ones  with  arbitrarily  large  gradients.  Of  course, 

waves  may  exhibit  both  properties. 

Numerical  computations  in  the  related  problem  of  the  interface  between  two  unbounded 

fluids  of  differing  densities  (cf.  [17],  [28])  show  that  there  are  branches  of  periodic 

waves  along  which  the  gradients  approach  infinity.  The  condition  (6.25)  necessitates 

w  being  large  if  the  vector  Vw  is  large.  However,  in  the  computations  cited  w 
x2  X1 

goes  to  infinity  and  thus  both  components  of  Vw  must  grow  unboundedly  along  the 
branch.  In  fact,  the  streamlines  steepen  to  the  point  of  having  vertical  tangents  and 
continue  past  this  configuration  to  one  of  "overhanging"  waves,  in  which  the  separating 
streamline  is  no  longer  a  graph  over  the  horizontal  axis.  Computations  are  underway  on 
the  solitary  wave  problem  studied  here  (cf.  [29]),  but  it  is  premature  to  predict  the 
large  norm  behavior.  We  will  show  that  if  waves  broaden  indefinitely,  maintaining 
uniformly  bounded  gradients,  then  a  solution  representing  an  "internal  bore"  with 
specified  amplitude  and  speed  (determined  by  (8.23)  and  (8.24))  must  exist.  If,  on 


the  other  hand,  there  were  e  sequence  e  S,  k  -  1,2,...  with  suplvw^l  ♦  •  as 

k  ♦  »,  one  could  expect  these  to  converge  to  a  limiting  wave  with  a  vertical  tangent  on 
the  interface. 


8. 1.  internal  Bores 

Let  (Xk,wk)  e  S,  k  -  1,2,  be  a  sequence  of  solutions  of  (2.36)-(2.39)  for  the 

density  p  9uppose  that  for  some  6  >  0,  Vw.  e  Q.  for  all  k  and  that  |wv|  .  ♦  m 

d  k  6  K1L2(T) 

as  k  ♦  ••  As  already  noted,  the  HSlder  estimates  (4.20)  and  (4.33)  for  the  gradients 

will  hold  uniformly  in  k.  Consequently  one  can  use  the  diagonal  method  to  select  a 

* 

subsequence,  still  denoted  ( X^ ,w^ )  for  which  X^  ♦  X  <  X^  and  w^  converges  in 
C0' '(B),  for  all  bounded  seta  B,  to  a  function  w(x^,x2)  e  C^'^tT)  n  C1,a(T*).  The 

a 

limiting  function  w  will  satisfy  all  of  the  estimates  derived  in  Lonmin  4.3  and  Theorem 

A 

4. 4 i  namely,  the  local  estimates  which  era  independent  of  R.  Naturally,  w  inherits  the 

A  A 

evenness  in  x 1  and  satisfies  3w/2x  1  <  0  for  x^  >  0.  Since  w  is  also  nonnegative, 

A  A 

lim  w(x1#x2)  -  (8.1) 

Xl"- 


exists.  To  examine  further  the  implications  of  this  convergence,  we  shift  it  to 

A  A  A  J 

B  “  [-1,1]  x  I  by  letting  w^fx^Xj)  "  wlXj  +  m,x2).  Now  w^  is  bounded  in  C  '  (B) 
uniformly  in  m,  according  to  the  HSlder  estimates  (4.20)  and  (4.33).  Since  w^ 
converges  in  C(B)  to  w<(  it  snist  also  converge  in  C®'  *(B  n  T*)  according  to  the 
Arsela-Ascoli  theorem  and  the  limit  w^  must  be  in  C  '  (I1)  as  a  function  of  x2 
alone. 

Now  let  8(x.|)  be  an  arbitrary  element  of  c"(-1, 1)  and  t,  an  arbitrary  element 

00  A 
of  C.(I) .  Using  the  product  8t  as  a  test  function  in  (4.39)  with  w  -  w  and 

U  01 

p  “  p and  letting  m  ♦  •  one  obtains 
•  a 


WV0' 


dw 

m 

dxT 


dx 

Ox, 


“  X(p_ 


P  +  )w-(0)T(0) 


(8.2) 
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which  is  the  week  form  of 

A 

d^pdf2(°'  •  (8-3) 

As  a  consequence  the  nonnegative,  continuous  function  w_  will  be  linear  where  ■=  0, 
that  is,  on  l“  and  on  I+. 

A  A 

Suppose  w^  =  0.  That  is,  suppose  that  w  -*■  0  as  x1  ♦  ■  .  Then  as  x^  ♦  »  the 
functions  w^  defined  at  the  outset  of  this  discussion  must  converge  to  zero  uniformly 

OB 

in  k  as  elements  of  L  (I)  depending  on  x^>  Their  gradients  must  converge  to  zero  in 
the  same  sense,  according  to  the  local  estimates  (4.20)  and  (4.33).  Suppose  X^  from  the 
pair  (\C'W)C)  converges  to  X  <  X^  as  k  ♦  •».  In  Lemma  5.2  an  estimate  was  made  of  the 
"tail”  of  wj  that  is,  the  restriction  of  w  to  Tm  »  (m,»)  x  I.  There  a  bound  for  w 

1  o,1- 

in  Hq (T)  was  used  to  show  that  the  C  '  (T^)  norm  of  w  went  to  zero  with  m, 

yielding  the  estimate  (5.6).  Since  the  functions  w^  under  discussion  have  their  norms 

in  C8,1(Tb)  converging  to  zero  with  increasing  m,  uniformly  in  k,  the  estimate  (5.6) 

can  be  seen  to  hold  for  w  »  wk  with  the  constant.  K  in  (5.6)  depending  only  on  X  and 

on  the  6  which  serves  for  all  k.  A  bound  for  w^  in  HQ([-m,m]  x  I)  follows  easily 

from  the  estimates  in  Section  4  and  thus  the  w^  are  uniformly  bounded  in  hJ(T), 

contradicting  the  initial  hypothesis  that  [w^j  1  ♦  ■  as  k  +  ». 

.  H0(T) 

In  the  case  that  Xfc  ♦  X  «  X^  as  k  ♦  ■  a  similar  situation  prevails  as  regards  the 

behavior  of  w  and  Vw  for  large  x^.  It  suffices  in  the  proofs  of  Theorems  5.6  and  5.7 

0,  1  - 

to  have  the  norm  of  w  in  C  (T^)  approach  zero  as  m  ♦  •  and  this  we  have  for 
w  “  w^,  uniformly  in  k.  A  consequence  of  those  theorems  is  that  w  *>  0  and  that  wfc 
converges  to  w  «  0  in  Hq(T),  again  contradicting  the  initial  assumption.  Consequently 
W"  and  w  are  both  positive.  He  summarize  what  has  been  proven  in 
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satisfy 


Iheorem  8.  1.  Suppose  that  f*jc'wjc)  f  S,  k  -  1,2,..., 

(a)  |Vw  |  <  K,  for  all  k 

L  (T) 

and 

(b)  X.  ♦  X  and  |w.  |  .  ♦  *  as  k  ♦  ®  . 

h’(T) 

Then  w^  converges  to  wtx^x^)  in  C(B)  n  C^B  n  T*)  for  each  bounded  set  B.  The 

A 

function  w  is  positive  on  T,  nonincreasing  for  x^  >  0,  real-analytic  in 
T  -  {x2  “  0},  and  satisfies 

lim  w(x,,«)  ®w  (8.4) 

x 

in  C(l)  n  C0' ’(I*5  where  w^  is  positive  on  I;  linear  on  1+  and  I“,  and  satisfies 
(8.2).  Both  w  and  w  satisfy  (2.36),  (2.37),  and  (2.39). 


Note  that  ^  is  a  solution  of  the  flow  equations  (just  enumerated)  which  is 

A 

independent  of  x.  Using  the  nomenclature  of  Benjamin  (9)  we  shall  call  wn  a  conjugate 

A 

flow.  According  to  Theorem  8. 1  the  function  ^  must  have  the  form 

A(x2  +  h). 


B(x, 


1  +  h), 


Xj  (  I 

x2  e  I* 


(8.5) 


for  some  value  of  A  where,  by  continuity  Ah  ®  B(-1  +  h).  Equation  8.2  gives 

A 

X  ■  A(A),  where 
A(A) 


7p — p)  IXfT  t  p_(  1 - - — ~  P+(  1 - - — 2^  1  ’ 

ip_  p+j  ■‘*n  (1  +  A)  (1  +  B) 


(8.6) 


Since  3w  /3x,  -  B  >  -1  on  I+,  we  have  0  <  A  <  (1  -  h)/h  in  (8.6).  The  amplitude 
A  and  the  speed 

c  -  cA  -  (g/A(A))1/2  (8.7) 

are  as  yet  unknown.  To  determine  A  we  show  that  as  a  consequence  of  Theorem  8. 1  there 
is  a  flow  in  the  strip  which  connects  a  trivial  flow  to  a  conjugate  flow.  That  is,  there 
is  a  solution  of  (2.36),  (2.37),  and  (2.39)  which  has  the  following  asymptotic  behavior. 


-1 
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As  x  it  approaches  a  flow  with  velocity  vector  (c^,  0 )  in  a  fluid  which  has 

density  p_  for  -h  <  y  <  0  and  density  p+  for  0<y<1-h.  As  x  +  +“  the 
conditions  approached  are 

( i  ■  *  A  cA,0)  for  -h  <  y  <  Ah 

velocity  “  ■  (8.8) 

((1  -  V^Ah)  CA'0)  for  Ah  <  y  <  1  -  h 
according  to  the  relations  (2.43). 

Theorem  8.2.  Suppose  the  hypotheses  of  Theorem  8.  1  are  satisfied.  Then  there  is  a 
solution  ( X,  j  (x^x^ ) )  of  (2.36),  (2.37),  and  (2.39)  with  the  following  properties: 

(a)  j  is  real  analytic  in  T  -  {x2  «■  0}  and  is  in  C1,a(T*)  n  C8'  1  (T )  for  some 
a  >  0. 

(b)  3  j/3x 1  >0  in  T. 

0.1- 

The  element  j(x1#«)  of  C  (I)  satisfies 

(c)  lim  |j(x  ,.)|  n  ,  _  “  0 

x,—  C0,1(I) 

and 

(d)  lim  |  j(x.,«)  -  wj  Q  -  0 

x^-f*  C°'\l) 

A 

where  w^  Is  a  conjugate  flow  corresponding  to  some  value  of  A  >  0  and  X  ”  A (A). 

Proof .  It  can  be  readily  verified  that 

AW  “  Ad  ~  '2(p_3-ep~)  *  *  0(a2)  (8‘9) 

for  A  near  zero,  where  Xd  and  e  are  defined  by  (2.40)  and  (2.41),  respectively. 

Since  e  >  0,  by  assumption,  A(A)  <  X^  for  small  positive  A.  Since  A(A)  ♦  +»  as 
A  ♦  (1  -  h)/h,  it  is  clear  that  A(A)  is  not  one  to  one  on  the  preimage  of  the  interval 
(0,Xd),  the  set  of  interest  to  us.  Let  the  conjugate  flow  w^  in  Theorem  8.1  correspond 
to  a  value  A  >  0  in  (8.5).  From  the  form  of  A ( A )  it  is  clear  that  there  is  a  smallest 


Implicit  in  the  fora  of  the  conjugate  flow  (8.5)  ia  the  conservation  of  mass  in  the 


lower  fluid,  connecting  a  flow  of  height  ( 1  +  A)h  and  speed  ~  at  x  *•  +■  with 


flow  at  -*>  of  height  h  and  speed  cA.  A  similar  situation  prevails  in  the  upper 


fluid.  This  same  conservation  can  be  obtained  starting  with 


div(pq)  “  0 


(8.  13) 


from  Section  2  (cf.  (2.9)— (2.  10))  since  streamlines  are  integral  curves  for  the  field  q. 


Another  conserved  quantity  can  be  derived  from  the  momentum  equation  (2.8).  It  can 


be  expressed  by  the  exactness  of  the  form 


1/2  . 

p  ud<i  +  p  dy 


(8.  14) 


assuming  now  that  we  are  working  with  a  smooth  density  p.  To  see  the  exactness  express 


(8.  14)  as 


p,/2u  -|i  dx  +  p  1/,2u  •|4-  dy  +  p  dy  "  -puv  dx  +  (pu2  +  p)dy  . 
dx  dy 


By  using  div  if  -  0  and  Vp  ♦  q  •  0  one  may  show  that 


(-PUV)  -  (pU  +  p)x  -  -ptU^U  +  UyV)  -  px 


using  the  horizontal  component  of  the  momentum  equation  (2.8). 


Now  41  and  y  are  constant  on  the  "walls"  at  y  ■  -h  and  y  »  1  -  h,  and  if 


(8.  14)  is  integrated  around  the  boundary  of  the  rectangle  (a,b)  x  (-h,  1  -  hj  there 


results 


/  (p1/2u  +  ■  /  <p1/2«  +  p)j  <*y 

-h  3y  'x«a  -h  3y  'x-b 


(8. 15) 


That  is,  each  integral  appearing  in  (8. 15)  is  independent  of  the  x  coordinate  (cf. 


[8]).  If  the  expression  (2.15)  for  the  total  head  N  is  used  to  replace  p  in  (8. 15) 


and  the  relation  (2.  14)  is  used,  the  x-invariant  quantity  in  (8.  15)  can  be  expressed  as 


/  lH  +  I  (*y  "  ’•’x5  "  gpyldy 
— h 


(8.  16) 


Given  the  regularity  and  the  convergence  properties  of  the  flows  considered  up  to  this 


V"  .■ W. • 


point  it  ia  easy  to  verify  that  the  expression  (8.  16)  will  be  invariant  for  the  flow 
associated  with  j  ( x  x2 )  *  By  equating  the  values  of  (8.16)  at  xj  «  -»  and  x 1  “  +~, 
one  obtains  another  condition  on  A  in  the  conjugate  flow  j+,  from  which  A  can  be 
determined. 

Normalize  the  pressure  so  that  it  is  0  at  xj  “  — x2  -  0.  Then  H  has  the 
values 

"*“TPtc2  (8*17) 

in  the  upper  and  lower  fluids/  respectively.  As  xj  ■  x  ♦  t",  ♦  0  for  the  flow 

2 

corresponding  to  j  and  the  term  p  is  merely  the  density  times  the  square  of  the 
horizontal  velocity.  At  x  ■  *•  the  expression  (8.16)  is 

2  2 
■j  p  cSi  +  rj  p  c^  +  gp_  +  -j  p+c2(  1  -  h)  +  -j  p+c2(  1  -  h)  -  gp+  — — 2  (8. 18) 

where  c  -  cA  in  (8.7).  The  value  of  (8.16)  at  x  -  ■*«•  is  computed  using  the  velocities 

and  heights  in  (8.8).  The  result  is 


IP-C2d  1  ♦  A 


c\  (A*h2  -  h2) 


-  <tP 


(1  -  h)2  -  A2h2 


+  1  p,c2( t  -  h  -  Ah,  +  1  P,  ^Z±L.  -  gp,  V [ '  ’  n,2  -  "  1  •  <*•”> 

Considering  the  expression  (8.19)  as  a  power  series  in  A  there  is  a  constant  term  which 
is,  of  course,  (8.18).  The  linear  term  in  A  vanishes  and  what  remains  is 


2  2  2 

1  2.  A  Ah  .1  2  1 

2  P-C^  TTa  -  *>-  —  +  2  P+°  — 


2  2 

_  +gp  AV 

Ah  .  + 


(8.20) 


(1 - 

'  1  -  h' 

which  must  be  zero.  If  an  A2  is  cancelled  and  g/c2  set  equal  to  X  the  result  is 


X 


P_ 

^h( 1  +  A)  + 


(8.21) 


If  this  is  equated  to  the  expression  (8.6)  and  terms  of  zeroth  order  in  A  are  eliminated 
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the  result  is 


2h( 1  +  A)2  2( 1  -  h  -  Ah)2 


(8.22) 


If  on  A  is  cancelled  and  the  result  solved,  one  finds 


/p_  ♦  /p+ 


(/p _  +  /p+)  /p.  /p  + 

hr  +  t=-f 


where  e  is  the  parameter  defined  in  (2.41). 

If  the  expression  for  A  in  (8.23)  is  used  in  (8.21),  the  result  is 

a  </p_  +  /p+) 

X  ■  _  • 

p_  -  Pa 


(8.24) 


This  last  expression  is  strictly  less  than  Ad  when  e  ?  0  and  is  equal  to  Ad  when 

e  -  0. 

Next  we  give  the 
8.2.  Proof  of  Lemma  6.4 

For  a  smooth  density  p  approximating  pd  and  a  solution  (A ,w)  of  (2.36)-(2.38) 
the  equation  (2.36)  can  be  written 


-p(x2,(ttvl  +  v  dr  rrnr  t}  "  pf2(7w)  *  'Xp‘w  • 

Xj  1  1  Xj  2 

For  each  x2,  w(x^,x2)  has  its  maximum  at  x^  "  0  and  hence  w^ ^( 0,x2)  “  0  and 
wx1x/°'x2^  *  °*  If  the  last  equation  is  restricted  to  the  line  {x^  -  0),  then 

-  p(x2)f2(0,wx  (0,x2))  <  -Ap'(x2)w(0,x2)  (8.25) 


results.  If  (8.25)  is  multiplied  by  g(x2)  “  w(0,x2)  and  the  result  integrated  over 


viRssSsssjsSs 

I,  it  yields 


/  pf,(0,g1>g'  <  -X  /  p'g2  •  (8.26) 

I  I 


Consider  now  s  pair  (X,w)  e  G.  .  There  is  a  sequence  (X_,w  )  e  F_  ,  D  with  X  ♦  X 

o#k  n  n  iifO ;  it  n 

and  w_  ♦  w  in  C0#  1 C  [—  1, 1]  *  1)  as  n  ♦  •  and  (8.26)  holds  with  p  ■  p  , 
n  n 

q  ■  q  -  w„(0,x.>),  and  X  “  X  .  If  F  “  F( p)  is  any  real  function  defined  for  p  >  -1 
n  n  *  n 

satisfying 

F(p)  <  f2(0,p)p  (8.27) 

for  -1+6<p<6-1,  then  it  follows  that 


X 


4 


inf 

g«Hg(i) 

-1+6<g'<4” 1 
0<g(0)<(1-h>(1-6) 


/  Pdy(9’) 
(p_  -  p+)g2(0) 


(8.28) 


is  a  lower  bound  for  X.  Note  that  g(0)  -  w(0,0)  >  0  for  any  nontrivial  element  of 
G.  ,  for  the  maximum  principle  applied  in  T+  and  in  T“  shows  w  has  its  maximum 

0  |K 

at  (0,0). 

The  function  F(p)  -  f(0,p)  is  easily  seen  to  satisfy  (8.27)  and  moreover,  is  convex 
for  p  >  -1,  having  a  second  derivative  equal  to  (1  ♦  p)“3.  Let  fj(p)  be  strictly 
convex  for  all  p  e  *»  coincide  with  f(0,p)  for  -1+6<p<6  ,  and  have  quadratic 

growth  near  p  «  *•».  One  can,  for  example,  extend  f °*P)  continuously  outside  of 
[-1  +  6,6  *)  to  be  constant  for  p  <  -1  +  6  and  for  p  >  6  '•  Two  integrations  will 
then  yield  a  suitable  f g.  The  infiinum  in  (8.28)  will  be  unaffected  by  replacing 
F(g* )  »  f (0,g* )  by  fg(g').  After  replacement  the  infimum  can  only  be  made  lower  by 
removing  the  condition  on  g*.  Thus 
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geH^(X) 

0<g( 0 )< ( 1— h  > ( 1—6 ) 


(p 


P+>g  (0) 


is  a  lower  bound  for  X. 

For  each  constant  y  >  0  the  set 


is  2m  affine  subspace  of 


Q'lge  H’(I)|g(0)  -  y> 
Hg(I)  and  the  convex  functional 
J(g)  -  /  Pdfs(g'> 


is  therefore  weakly  lower  semi-continuous  on  Q.  Since  J  is  nonnegative,  there  is  a 
g0  e  Q  at  which  J  assumes  its  minimum.  For  n  <  cj(l  ),  g,.,  +  tn  €  Q  and  a  standard 
variational  argument  shows 

/  Pd^O^'  *  0  * 


From  the  classical  theorem  of  IXiBois-Reymond  it  follows  that  g£  is  constant  on  I-- 
Likewise  it  is  constant  on  I*.  Thus  gg  has  the  form  (8.5)  of  a  conjugate  flow  with 
g0(0)  -  Ah  <  (1  -  h)(1  -  6).  For  such  a  function  the  quotient  in  (8.29)  takes  the  form 


1 

2 


P 


A2  1  A2h2 

TTa  h  +  2  P+  1 - 


h  -  Ah 


(P+  -  p_)A2h2 


(8.30) 


which  is  precisely  half  the  expression  in  (8.21).  At  a  minimum  A  satisfies  (8.22)  and 
the  corresponding  minimum  value  is  half  of  that  in  (8.24),  giving  the  bound  (6.12). 
q.e.d. 


A  better  lower  bound  for  X  could  be  found.  The  function 


d  (2  +  p) 

f2'°'P>P  “  -  2 

2d  +  Pr 


(8.31) 
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*5 


is  convex  on  -1  +  6<p<2  for  each  4  >  0  and  can  be  extended  outside  that  range/  ae 
before,  to  a  globally  defined  convex  function  (p)  which  grows  quadratically  at  «• • 

One  can  now  use  an  argument  like  that  just  given*  The  expression  in  (8.31)  can  be  used 
for  P(p)  in  (8.28)  with  the  restriction  -1  +  4  <  g'  <2.  In  (8.29)  one  replaces  f, 
by  and  finds  a  piece-wise  linear  minimizer  g  having  a  derivative 


g*(x2) 


1  -  h' 


x2  <  0 
x2  >  0 


(8.32) 


Note  that  K  e  (0,  --  ~  ~)  since  g*  >  -1.  If  h  >  1/3,  so  that  (1  -  h)/h  <  2,  then 
g*  lies  in  the  range  for  which  the  Integrand  in  (8.29)  is  the  expression  in  (8.31). 


Corresponding  to  (8.30)  one  obtains  an  expression  to  minimize  for  A  e  (0,  —  ) 


Unfortunately,  the  minimizing  value  of  h  satisfies  a  fourth  order  equation.  Without 
giving  the  details  we  merely  note  that  if  p+  ■  2,  p_  -  1,  and  h  -  1/2,  so  that 
(1  -  h)/h  <  2,  the  following  numerical  values  are  obtained.  The  bifurcation  point  is 

“  6  and  Lemma  6.4  gives  the  lower  bound  2.91,  whence  2.91  <  \  <  6.  The  minimization 
using  (8.31)  yields  the  better  estimates  5.81  <  A  <  6. 
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